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THERMIONIC  DIODE  CURRENTS  IN  TRAIBVERSE  MAGNETIC  FIELDS* 
R.  J.  Zollveg  and  M.  Gottlieb 


Introduction 

The  transport  of  electrons  from  cathode  to  anode  of  a  tbenalonic 

diode  is  Impeded  by  the  presence  of  a  transverse  magnetic  field.  If  the 

electron  flow  is  in  a  vacuum  the  well-known  magnetron  effect  takes  place. 

In  a  gas  the  electron  diffusion  and  mobility  are  not  reduced  as  greatly 

by  the  transverse  magnetic  field  because  the  scattering  of  electrons  by 

gas  atoms  randomizes  the  direction  of  electron  motion  with  each  collision. 

For  thermionic  power  converter  diodes  it  has  been  recognized 

that  the  magnetic  field  induced  by  the  diode  current  itself  can  be 

sufficient  to  seriously  reduce  cell  current  thus  producing  a  practical 

upper  limit  to  diode  current-carrying  capacity.  Schock^  has  made 

extensive  calculations  on  the  effects  of  self-induced  magnetic  fields 

2 

for  thermionic  diode  ciirrents  in  a  vacuum.  Houston  has  made  a  few 
measurements  in  cesium  vapor.  More  extensive  and  precise  data  are  needed 
for. converter  design  calculations.  When  this  work  was  undertaken  it  also 
appeared  that  it  might  be  possible  to  obtain  Infomatlon  about  the 
collision  frequency  and  elastic  scattering  cross  section  of  electrons  in 

3 

neutral  cesium  since  on  the  basis  of  relations  derived  by  Langmuir  and 
Coblne  ,  the  measurement  of  cell  current  in  a  gas  should  give  a  measure 


Ihls  work  was  sponsored  by  Aeronautical  Systems  Division,  Air  Force 
Systems  Command,  United  States  Air  Force,  in  part. 

1.  A.  Schock,  J.  Appl.  Phys.  ^  1978  (i960). 

2.  V.  C.  Wilson  ARS  Space  Power  Systems  Conf.,  Santa  Monica,  California, 
September  27-30,  i960. 

3.  I.  Langmuir,  Phys.  Rev.  ^  I656  (1931). 

4.  J.  D.  Cobine,  Gaseous  Conductors.  McGraw-Hill  Book  Co.,  Inc. 

New  York  (1941). 
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of  the  mobility.  Our  results,  however,  show  that  although  the  calculated 
cross  sections  are  of  the  right  order  of  magnitude,  difficulties  most 
likely  associated  with  back  diffusion  of  electrons  to  the  emitter  prevent 
the  precise  determination  of  scattering  cross  sections.  Nonetheless  the 
data  should  be  of  considerable  value  in  connection  with  converter  design 
since  the  combined  effects  of  electron  transport  and  back  diffusion  In 
a  transverse  magnetic  field  In  cesium  vapor  eire  not  at  present  calculable 
from  first  principles. 


follows : 


An  empirical  result  useful  for  converter  design  is  obtained  as 


,  ,  k.a  X  10‘5 


where  is  the  current  in  the  absence  of  the  magnetic  field,  I  with 
magnetic  field,  d  Is  emitter-collector  spacing  in  Inches,  B  is  magnetic 
field  in  gauss  and  p  Is  Cs  pressure  In  mm  Eg. 

II.  Experimental  Technique  and  Results 

Cell  current  measurments  were  made  with  cells  of  two  different 
structures  as  shown  in  Figure  1.  Both  cells  utilized  a  O.CX}^  inch  by 
O.OSO  inch  by  1  Inch  tungsten  ribbcn  filament  heated  directly  with 
halfwave  rectified  A.C.  power.  The  anode  in  one  case  consisted  of  plane 
peuredlel  copper  plates  with  cathode -euaode  spacings  of  0.020  inches.  The 
anode  in  the  other  case  was  a  0.^00  inch  diameter  stainless  steel  cylinder 
which  thus  had  an  average  cathode-anode  spacing  of  0.225  Inches,  about 
an  order  of  magnitude  larger  than  the  first  cell.  Ihe  tubes  were  mounted 
in  an  oven  inside  a  solenoid  magnet  the  field  of  which  could  be  continually 
varied  up  to  about  350  gauss. 

The  circuit  shown  in  Figure  2  was  used  for  the  measurements. 

Cell  currents  were  determined  by  the  voltage  drop  across  a  known  resistance 
during  the  time  that  the  filament  voltage  was  zero,  using  a  calibrated 
oscilloscope.  To  reduce  the  scatter  of  the  data  associated  with  changes  In 
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filament  emission  a  switch  was  used  that  simultemeously  closed  the 
magnet  current  circuit  and  replaced  the  cell  load  resistance  by  8Ui 
alternative  resistsince  box.  Thus  the  load  resistance  when  the  magnetic  field 
was  present  could  be  increased  to  compensate  for  the  reduced  cell  current. 

In  this  way  the  cell  voltage  and  internal  fields  remained  constant.  A 
Clare  Hg-wetted  relay  synchronously  opened  the  anode  lead  during  the 
filament  heating  pulse.  A  battery  was  added  to  the  circuit  to  supply 
positive  potentials  to  the  anode.  The  required  cesium  vapor  pressure 
was  established  by  emersing  a  side  tube  containing  the  cesium  ampule 
in  a  Sn-Pb  eutectic  bath  which  could  be  controlled  to  within  about  1°G. 

The  rest  of  the  envelope  was  heated  to  a  hi^er  temperature  to  prevent 
Cs  condensation.  The  communicating  apertures  between  the  diode  and 
the  tube  enve‘'ope  were  sufficient  that  estimates  of  the  themELL  trans¬ 
piration  effect  lead  to  the  conclusion  that  appreciable  pressure 
differences  do  not  exist  between  the  cathode -anode  region  and  the  cesium 
reservoir  over  most  of  the  region  of  interest. 

The  results  of  the  measurements  were  found  to  depend  upon 
filament  temperature  and  cell  potential  as  well  as  expected  dependence 
upon  cesium  pressure  and  magnetic  field  strength.  With  the  larger 
spacing  cell  the  results  with  lower  filament  temperatures  and  higher 
filament  temperatures  were  not  far  different  while  intermediate  tempera¬ 
tures  gave  a  greater  reduction  in  cell  current  with  magnetic  field. 

Since  at  the  higher  temperature  used  one  would  expect  sufficient  ions 
to  be  formed  by  contact  ionization  at  the  cathode  to  neutralize  the 
electron  space  charge  and  since  the  space  charge  is  not  as  severe  a 
problem  at  the  lowest  cathode  temperatures  where  electron  emission  is 
smaller,  the  higher  current  reductions  for  intermediate  temperatures 
appear  to  be  associated  with  additional  space  charge  effects  associated 
with  the  reduced  electron  mobility  in  a  magnetic  field.  In  Figure  5 
are  plotted  some  of  the  results  for  the  larger  spacing  cell.  We  have 
plotted  data  only  for  the  higher  emitter  temperature  (2000°k)  where  the 
effects  of  space  charge  should  be  negligible.  It  should  be  remembered, 
however,  that  the  magnetic  field  greatly  decreases  the  electron  mobility 
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without  a  comparable  decrease  in  the  cesium  ion  mobility.  The  ratio 
of  the  current  in  the  absence  of  the  magnetic  field  to  that  with  the 
magnetic  field  has  been  plotted  vs.  the  square  of  the  magnetic  field 
strength  for  three  cesixam  reservoir  temperatures  and  two  cell  potentials. 
Note  that  the  current  reductions  are  somewhat  less  for  higher  anode 
potentials.  Included  is  the  curve  for  the  cesium  reservoir  at  room 
temperature --where  electron-cesiian  atom  scattering  is  negligible. 

Because  of  the  contact  potential  difference,  the  curve  obtained  with  a 
battery  potential  of  I.56  volts  probably  corresponds  to  approximately 
a  zero  field  condition  between  the  electrodes. 

Figure  4  shows  the  results  for  the  closer  spacing  cell.  All 
of  the  1q/1  vs.  B  curves  in  Figure  4  are  linear,  except  the  one 
corresponding  to  =  250°C  with  a  filament  temperature  of  609°K  emd 
+8.8  volts  cell  potential.  Most  of  the  results  with  this  cell  were 
Independent  of  cell  potential  as  shown  in  Figure  ^  until  potentials  suf¬ 
ficient  to  cause  inelastic  scattering  and  perhaps  volume  ionization  are 
reached.  It  is  evident  that  there  is  considerable  variation  in  the 
slopes.  The  decrease  in  slope  with  Increase  in  cesium  pressure  is  less 
than  might  be  e^qiected  as  discussed  in  the  next  section. 

In  the  low  magnetic  field  region  where  qualitative  con3>arlson 

is  possible  the  results  are  in  reasonable  agreement  with  those  of 
2 

Houston  . 

III.  Interpretation 

The  calculations  of  Schock^  for  electron  treuasport  perpendicular 
to  the  magnetic  field  in  a  vacuum  grossly  overestimate  the  effect  of  the 
magnetic  field  when  cesium  vapor  Is  present.  One  expects  the  electron 
transport  in  this  case  to  be  dominated  by  the  relations  for  the  mobility 
and  diffusion  in  the  magnetic  field. ^ 


_ 2 _  D  =  _ 2 _ 

(l+flcpT^)  ’  (l+tU^T^) 


(1) 


where  and  are  the  values  in  the  absence  of  the  magnetic  field, 

o)  =  —  B  and  t  is  the  period  between  collisions.  Equations  (l)  are 
m 


5.  W.  P.  Allis,  Handbuch  der  Physik  XXI,  585,  Springer-Verlag,  Berlin  (1956). 
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strictly  valid  only  if  the  collision  frequency  is  independent  of 
electron  velocity  hut  are  also  a  good  approximation  in  other  cases  for 
sufficiently  high  magnetic  fields. 

In  cases  where  there  is  no  space  charge  and  only  elastic 
collisions  Langmuir^  found  that 


where  is  the  saturation  current,  \  is  the  electron  mfp,  x  is  plate 

separation,  is  cell  potential  and  is  the  initisd  electron  energy. 
Since  both  the  diffusion  constant  and  mobility  are  proportional  to  the 
mean  free  path  one  would  e:q>ect  in  this  case  that 


^0 

r 


1  + 


(2) 


where  1^,  are  the  values  without  anc\  I,  in  the  presence  of  a 

transverse  magnetic  field.  The  question  of  space  charge  limited  currents 

4 

in  the  presence  of  a  scattering  medium  has  been  considered  by  Coblne 

7 

end  in  greater  detail  by  Shockley  end  Prim  .  Where  initial  velocities 
can  be  neglected  they  obtain 

9K 

I  =  - jr — 

8  x^  (3) 


for  the  plane  parallel  case  with  electrode  spacing  x,  a  result  analogous 
to  the  Lengmuir-Child' s  Law  for  the  vacuum  diode.  If  Eq.  (3)  were  to 


7.  W,  Shockley  end  R.  C.  Prim,  Phys.  Rev.  22.^  T53,  (l955). 
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apply  then  Eq.  (2)  would  also  be  valid.  Eq.  (5)  does  assume  that  the 

electron  mobility  is  independent  of  energy,  as  has  recently  been  pointed 

8 

out  be  Formsui°  who  derives  another  expression  valid  for  rare  gases  under 
certain  conditions. 

Our  measurements  shown  in  Figures  2  and  5  are  approximately 

I  2 

consistent  with  a  relation  of  the  form  _o  =  1  +  K  B  as  expected 

from  Eq.  (2).  But  there  is  some  scatter  in  the  data  and  some  depetrture 

I  2 

from  linearity.  Furthermore,  the  slope  K  of  the  curve  _o  vs  B  ,  which 

2  II 

should  be  proportional  to  t  and  hence  to  is  found  not  to 

vary  in  the  proper  way  with  pressure.  There  is  also  a  dependence  on 
electrode  spacing.  The  most  likely  source  of  this  discrepancy  would 
appear  to  be  the  back  diffusion  of  electrons  from  the  plasma  to  the 
electrodes  which  might  reasonably  be  expected  to  depend  upon  pressure, 
magnetic  field  and  spacing.  While,  in  principle,  measurements  of  this 
sort  might  be  used  with  an  assumed  mobility  dependence  to  determine  the 
back  diffusion  contribution,  it  is  felt  that  the  results  are  not  siaffi- 
ciently  precise  to  merit  this. 

For  converter  operation  the  results,  of  course,  are  directly 
applicable.  The  results  in  Figures  3  and  4  should  give  reasonable 
estimates  within  perhaps  a  factor  of  2  of  the  magnitude  of  the  effect. 
For  other  spacings  and  pressures,  estimates  can  be  made  by  use  of  the 
empirical  result 


d^*^  B^ 

where  d  is  the  emitter-collector  spacing  in  inches,  B  is  in  gauss  and 
p  is  the  pressure  in  mm  Hg.  The  reduction  in  cell  current  should  be 
diminished  somewhat  at  higher  ion  concentrations  because  of  the  larger 
scattering  cross  section  of  Cs^  compared  to  neutral  Cs. 


»  1  +  X  10'^ 

I 


8.  R,  Forman,  Phys.  Rev.  123,  1557  (l96l). 
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A  summary  of  the  available  data  on  electron  scattering  by  Cs 
atoms  is  shown  in  Figure  6.  The  solid  curve  shows  Erode 's  data.  The 
circles  show  the  results  as  originally  given  by  Boeckner  and  Mohler,  but 
in  a  later  paper  they  indicated  that  their  electron  density  measurements 
were  in  error  at  the  higher  pressures  and  the  squares  show  these  results 

g 

corrected  by  Phelps  using  the  results  given  by  Mohler 's  paper.  Steinberg's 
data  was  also  used  to  calculate  collision  frequencies  and  his  values  at 
the  two  lower  Cs  pressures  are  shown  here  although  his  electron  densities 
presumably  also  may  be  in  error.  It  should  be  noted  that  the  data  of 
Erode  need  not  necessarily  agree  with  that  of  Mohler  since  Itohler's  results 
are  for  a  mobility  measurement  while  Erode 's  measurements  would  be  sensitive 
to  very  small  angle  scattering.  One  might  eaqpect  Erode 's  results  to  yield 
a  larger  collision  frequency  than  Mohler 's  contrary  to  Figure  6. 

We  wish  to  thank  A.  V.  Phelps  for  several  discussions. 


9-  A.  V.  Phelps,  Private  communication. 
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Fig.  2-Circuit  for  measuring  dependence  of  cell  current  upon  magnetic 

field  strength 
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Fig.  6 


Available  data  on  electron  scattering  in  cesiiom  vapor. 
Collision  frequency  Cs  density  per  cm3  is  plotted  vs. 
electron  energy. 
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FABRICATION  OF  ELECTRON  BEAM  HEATED  BUTTON-CATHODE  CELL 


R.  J.  Zollweg  -  M.  Gottlieb  C.  R.  Taylor 


Fabrication  has  recently  been  brought  to  completion  on  a 
cell  in  which  the  cathodes  are  indirectly  heated  by  electron  beam 
bombardment.  The  assembled  cell  is  shown  in  the  photograph  in  Figure  1. 
It  is  an  all  metal  system  which  is  demountable  at  two  shear  gasket  seals 
into  three  sections.  A  photograph  of  the  disassembled  cell  components 
is  shown  in  Figure  2.  The  two  pieces  on  the  right  are  joined  together 
in  the  final  braze  and  form  the  collector  assembly.  The  section  on  the 
left  is  the  electron  gun  assembly.  It  will  be  continuously  pumped  by 
an  auxiliary  vacuum  system.  The  electron  gun  filament  is  mounted  on 
two  of  the  high  voltage  lead  throughs,  and  the  other  two  may  be  used 
for  a  focusing  cone  for  the  beam,  if  necessary.  The  collector  is  the 
bottom  surface  of  the  central  cylinder  of  the  right  hand  structure. 

It  is  made  of  stainless  steel,  but  plates  of  other  anode  materials  may 
be  fastened  over  it.  The  collector  cylinder  is  hollow,  and  a  heater 
may  be  inserted  inside  to  regulate  the  temperature  of  the  anode.  The 
position  of  the  anode  can  be  adjusted  by  means  of  the  flexible  bellows, 
which  is  driven  by  the  three  micrometer  screws.  The  cathode  and  anode 
are  insulated  from  each  other  by  the  alumina  seal  above  the  bellows, 
and  by  the  ceramic  posts  attached  to  the  three  micrometer  screws.  The 
cathode -anode  area  is  viewed  through  a  sapphire  window  just  beneath 
the  bellows.  The  sapphire  disc  is  brazed  to  a  kovar  cylinder  with  BT 
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solder  using  titanium  hydride  as  flux.  The  cathode  temperature  and 
electrode  spacing  are  measured  through  this  window.  The  cathodes  to 
he  used  in  this  cell  are  in  the  form  of  buttons ^  l/2"  in  diameter 
and  l/8"  thick.  These  buttons  are  brazed  onto  a  0.010"  thick  tantalum 
foil.  A  photograph  of  the  entire  cathode  assembly  is  shown  in  Figure  3* 
The  cathode  in  the  photograph  is  tungsten,  and  the  braze  was  accomplished 
by  melting  a  niobium  foil,  whose  melting  temperature  is  lower  than 
that  of  either  tungsten  or  tantalum,  between  the  button  and  the  tantalum. 
The  niobium  alloys  with  both  metals,  and  forms  a  good  thermally 
conducting  bond  between  them.  Other  brazing  techniques  may  have  to  be 
used  with  other  cathode  materials.  The  sealing  gasket  is  the  outer, 
flat  edge  of  the  tantalum  support  foil,  which  is  sufficiently  ductile 
to  form  a  shear  seal.  Each  such  seal  can  then  be  made  only  once  with 
one  support  foil.  The  cell  will  be  connected  to  a  vacuum  system  through 
a  Granville -Phillips  brakeable  valve  for  processing. 


i 
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THE  LIFETIME  AND  EFFICIENCY  OF  A  THERMIONIC  ENERGY  CONVERTER 
L.  S.  Rlcbard8on>  A.  E.  Fein,  M.  Gottlieb,  G.  A.  Keneny  and  R.  J.  Zollveg 

I.  Introduction 

The  design,  construction,  and  operation  of  all  devices  are 
dictated  by  their  required  lifetlaes.  !I%ils  requlresient  varies  from 
minutes  for  a  rocket  motor  to  yesurs  for  a  coonerclal  power  plant. 

The  mechanisms  of  failure  must  be  considered  In  order  to  predict  this 
lifetime.  In  a  thermionic  energy  converter,  the  basic  cause  of  failure 
Is  the  evaporative  transfer  of  cathode  materiads.  Other  causes  such  as 
corrosion,  mechanical  failure,  or  radiation  damage  can  In  most  eases 
be  prevented  by  ingenious  design.  Evaporative  loss,  however,  is  an 
Inherent  limitation  because  of  the  hi^  operating  temperatures  required. 

To  predict  the  lifetime  of  a  thermionic  energy  converter,  only  the 
allowable  loss  and  the  rate  of  evaporative  transfer  must  be  known.  The 
adlowable  loss  Is  determined  by  the  exact  smchanlsm  of  failure.  Whisker 
formation,  changes  In  the  thermal  emisslvlties,  auid  chatnges  in  the  work 
functions  sure  examples  of  such  isechsualsms,  any  of  which  may  be  caused  by 
evaporation  emd/or  deposition.  The  rate  of  evaporative  transfer  Is 
dependent  upon  the  geometrlced  configuration,  the  operating  temperatures, 
the  materials  of  construction,  smd  the  Internal  pressure  of  the  device. 

An  energy  converter  must  possesB  a  useful  energy  output  as 
well  as  a  useful  lifetime.  In  order  to  compeure  the  operation  of  different 
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energy  converters,  some  criterion  of  comparison  must  be  chosen.  Some 
possible  criteria  are  the  efficiency,  the  output  power,  or  the  power 
to  weight  ratio.  Each  of  these  is  important  in  certain  applications; 
for  example,  a  stationary  land-based  power  station  is  usually  an 
efficient  low  cost  device  while  a  space  power  plant  is  usually  a  ll£^t 
device.  In  this  paper,  the  efficiency  has  been  taken  as  the  criterion. 

Iliis  choice  was  made  for  several  reasons:  l)  The  efficiency  is  eQ.ways 
of  interest  to  the  designer  regeordless  of  the  application  of  the  device. 

2)  Calculation  of  the  power-to-weight  ratio  requires  a  detailed  knowledge 
of  the  heat  source  and  sink,  and  is  thus  dependent  upon  considerations 
beyond  the  scope  of  this  paper.  This  point  will  be  discussed  further  in 
a  later  section.  3)  The  output  power  and  other  possible  criteria  are 
more  specifically  related  to  particular  applications  than  is  the  efficiency. 

A  number  of  authors  have  calculated  efficiencies  of  thermionic 

energy  converters  using  a  number  of  different  approximations. While  the 

efficiency  is  a  function  of  a  leurge  number  of  variables,  in  principle  most 

of  these  variables  are  determined  by  a  relatively  small  number  of  parameters. 

If  the  materials  of  construction  and  their  properties  are  known,  then  the 

efficiency  is  a  function  only  of  the  cathode -to-anode  spacing,  the  cathode 

and  anode  temperatures,  the  means  of  space  charge  neutralisation,  and  the 

lead  and  load  resistances.  In  most  converters  presently  under  investigation, 

cesium  vapor  is  used  for  space  charge  neutreilizatlon,  and  frequently  also 

for  modification  of  the  cathode  and  anode  work  functions,  as  recently 

(2) 

patented  by  G.  R.  Feaster.'  '  In  this  case,  the  efficiency  is  also  a 
function  of  the  cesium  pressure.  In  all  but  the  most  recent  calculations 
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of  efficiency  referred  to  above,  the  interrelationships  of  work  functions, 
emissivities,  temperatures,  etc.,  have  been  largely  Ignored.  Furthermore, 
few  of  the  csdculations  have  considered  variations  in  load  resistance  or 
the  presence  of  back  current.  None  of  the  calculations  have  atten^ted 
to  optimize  the  values  of  anode  and  cesium  temperature  while  simultaneously 
optimizing  lead  and  load  resistance.  A  recent  report  (3)  summarized  a 
method  for  this  simultaneous  optimization. 

Both  lifetime  and  efficiency  can  be  shown  to  be  functions  of 
the  seune  design  parameters:  operating  temperatures,  internal  pressure, 
materials  of  construction,  and  geometrical  configurations.  For  a  converter 
built  of  a  given  material,  and  operating  in  the  most  efficient  manner  possible, 
lifetime  can  be  shown  to  decrease  with  Increasing  efficiency.  The  bedanclng 
of  the  desires  for  high  efficiency  and  long  lifetime  ret^uires  the  knowledge 
of  how  these  quantities  change  as  the  design  parameters  are  changed.  In 
this  paper,  a  proposed  method  for  calculating  these  quantities  is  discussed. 
The  relationship  of  efficiency  and  lifetime  derivable  from  these  ccd.culatlons 
is  also  presented. 

The  methods  by  which  lifetime  and  efficiency  can  be  calculated 
are  described  in  Sections  II  and  III  respectively.  The  relationship 
between  these  two  properties  is  discussed  in  Section  IV. 

II.  Lifetime 

A  simple  approach  is  to  estimate  lifetime  by  dividing  the  allowable 
surface  recession  by  the  rate  of  evaporation  into  vacuum.  Since  no  measure¬ 
ments  of  the  allowable  surface  recession  have  appeared  in  the  literature,  a 
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value  of  2.5  microns  was  used  as  a  conservative  estimate.  This  value, 

4 

equivalent  to  about  10  layers  of  atoms,  is  small  compared  to  the  spacing 
of  most  converters  presently  under  consideration.  In  a  later  section  the 
exact  value  used  will  be  shown  to  be  unimportant.  I^igure  1  shows  the 
results  of  this  calculation  for  molybdenum,  columbium,  tantel-um,  and 
tungsten  (4). 

Several  important  effects  occur,  particularly  in  the  cesium 
diode,  that  are  neglected  in  this  calculation.  !nie  presence  of  a  partial 
cesium  coverage  may  Increase  or  decrease  the  rate  of  loss.  The  presence 
of  contaminant  gases  (such  as  water  vapor)  in  the  device  may  also  change 
the  rate  of  loss.  Experimental  measurements  under  very  carefully  controlled 
conditions  are  necessary  to  determine  the  magnitudes  of  these  effects. 
Another  effect  of  considerable  importance  is  the  reduction  of  evaporative 
transfer  caused  by  back- scattering  from  the  cesium  atmosphere  present  in 
the  anode  to  cathode  gap.  The  magnitude  of  this  effect  can  be  estimated 
by  a  relatively  simple  method. 

A  first  approximation  is  to  assume  that  all  atoms  that  strike 
the  cold  surface  (anode)  stick  and  that  the  partial  coverage  of  cesium  has 
no  effect  on  the  evaporation  of  atoms  from  the  cathode. 

For  relatively  large  spaclngs,  the  rate  of  loss  is  limited  by 
the  rate  of  diffusion  through  the  cathode -to-anode  gap.  If  the  diffusion 
constant  D  does  not  change  within  the  gap,  then  the  flvK  R  of  metal 
transported  is 


R  =  D 


Af 

d 


(1) 
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where  ^  f  is  the  difference  in  density  of  evaporated  atoms  at  the  cathode 
and  anode,  and  d  is  the  width  of  the  cathode-to-anode  gap.  Ihe  evaporative 
flux  at  any  temperature  T  is  given  by 


\  =  PVA  (2) 

where  is  the  density  and  V  is. the  arithmetic  average  velocity  of  the 

evaporating  species. ^^^Slnce  the  pressure  and  tenperature  of  evaporated  atoms 

adjacent  to  the  electrodes  are  nearly  equal  to  those  characteristic  of  the 

electrodes,  equsd  to  the  difference  of  densities  characteristic  of 

cathode  and  anode.  Since  the  cathode  temperature  is  considerably  higher 

than  the  anode  temperature,  the  density  at  the  cathode  ^  is  much  greater 

c 

than  that  at  the  anode  P Thus  is  equal  to  and  the  transport 
flux  is 


p  4D  ^vc 
^  “  V  d 


(5) 


where  R^^  is  the  normal  evaporative  flux  at  the  cathode  temperature. 

The  quantity  D/V  has  the  dimensions  of  length  and  is  proportions^, 
to  the  mean  free  path,  Hard  sphere  calculations  (6)  show  that  D/V  equals 
\/5. 

Thus, 


R  »  (V3)(V<i)  Rye  (^) 
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The  evaporative  transport  problem  is  analogous  to  the  neutron 
scattering  problem  in  a  non-absorbing  medium,  where  is  aneilogous  to 
the  incoming  neutron  intensity  and  R  to  the  intensity  transmitted.  C.  W. 
Maynard  has  solved  this  neutron  scattering  problem  numerically  for  small 
spacings  (7).  Examination  of  his  results  between  1/4  and  4  mean  free 
paths  shows  that  the  equation 


R  -  R^^(l  +  5d/4\) 


(5) 


fits  the  results  very  accurately.  For  large  spacing  this  expression 
approaches  that  of  equation  4. 

A  more  general  approach  js  to  consider  the  space  as  a  succession 
of  reflecting  and  transmitting  slabs,  llie  transmission  through  each  of 
these  slabs  must  have  the  same  functioned  relationship  to  thickness  as 
the  transmission  through  the  assembly.  If  the  multiple  reflections  are 
summed  up,  assuming  isotropic  behavior,  the  functioned  relationship  is 


f(I ».)  = 

i  ^ 


Z( 


l-f(x^) 


-1 


'(6) 


where  is  the  thickness  of  the  ith  slab  and  f(Xj^)  is  the  fraction 
transmitted  through  a  thickness  x^. 

It  can  easily  be  shown  that  the  only  function  which  obeys  this 
relationship  is 


-1 


f(x)  =  (1  +  ax) 


(7) 
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This  implies  that  R  ■  — is  the  general  form  of  the  solution  to  the 
transport  equation.  For  a  sufficiently  large  d,  the  constant  term  can  he 
neglected  in  comparison  vlth  ad,  thus  reducing  the  equation  to  the  diffusion 
case  where  a  is  3/^\.  The  scattering  approach  of  Maynard  slLso  led  to  the 
same  value  for  the  constant.  Both  approaches  axe  based  on  all  collisions 
being  of  the  hard  sphere  variety.  Thus  equation  3  is  the  genereil  solution 
for  transport  where  the  factor  3/4  will  be  different  if  the  collisions  do 
not  obey  the  hard  sphere  model. 

An  experlnental  verification  of  equation  3  was  attempted.  A  diode 

was  constructed  using  a  molybdenum  cathode  and  a  nickel  anode.  A  schr;matic 

diagram  of  this  cell  is  shown  in  Figure  2.  Four  different  anode  to  cathode 

spaclngs  were  used.  The  cathode  was  held  at  a  fixed  temperature  for  a 

measured  length  of  time,  after  which  the  nickel  anodes  were  analyzed  for 

total  molybdenum  content.  The  speetrographic  analysis  had  a  sensitivity 

of  about  2  micrograms  of  molybdenum.  Ihe  reproducibility  of  analysis  was 

estimated  to  be  about  13^  of  the  analyzed  value.  Table  I  summarizes  the 

results  of  the  first  test  of  this  cell.  The  ratio  of  the  fluxes  R  /R 

vc' 

is  clearly  the  same  as  the  ratio  of  the  eunounts, 

If  these  data  are  plotted  as  a  function  of  spacing  (Figure  3) 
equation  3  can  be  used  to  calculate  values  of  the  mean  free  path.  The  mecui 
free  path  deduced  is  about  8.6  microns  at  a  cesium  pressure  of  1  mm  of  Bg. 

No  other  measurement  of  mean  free  path  under  these  non- isotherms I 
conditions  was  found  in  the  literature.  TheoreticeO.  estimates  assuming 
isothemal  conditions  were  0.9  to  20  times  the  experimented  value,  depending 
on  the  assumptions  made. 
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since  the  mean  free  path  is  inversely  proportional  to  the  pressure 
and  is  not  a  sensitive  function  of  the  temperatures,  the  following  relationship 
was  assumed  to  hold  true  for  the  transport  flux 

R  -  (1+0.087  pd)‘^  (8) 

for  the  pressure  p  in  mm  of  Bg  and  the  spacing  d  In  microns.  The  value  of 
0.087  Is  based  upon  only  one  series  of  experiments  and  should  be  accepted 
only  with  considerable  reservation. 

From  equation  8,  It  is  possible  to  determine  lifetime  as  a  function 
of  cathode  temperature,  spacing,  axid  cesium  bath  temperature.  relative 
lifetime,  defined  as  the  ratio  of  the  lifetime  in  a  gas  filled  diode  to  the 
lifetime  In  vacuum.  Is  equal  to  In  Figure  4  the  relative  lifetime 

Is  plotted  as  a  function  of  spaclxig  and  cesium  bath  temperature  which  deter¬ 
mines  the  pressure.  Better  determinations  of  the  mean  free  path  as  a  function 
of  pressure  may  change  the  values  shown  In  this  figure,  but  the  general 
validity  of  the  figure  should  not  be  affected.  The  lifetime  will  obviously 
Increase  with  increasing  cesium  pressure  and  spacing. 
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Tteible  I  -  Experimental  Determination  of  Molybdenum  Transport 


Spacing  d  (ml crons) 
155 
226 
450 
450 


Amount  Transported  G(tig) 

19 

12 

6.2 

2.6 


Temperature  m  2100* K 

2 

Time  ■  49^200  seconds  Area  of  Cell  >  0.2  cm 
Evaporation  in  Vacuum,  0^  ■  250  pg  (Reference  4) 
Cesium  Temperature,  552* K 
Cesium  Pressure,  1.0  mm  cf 


-30 


III.  Efficiency 

The  basic  equation  used  for  the  calculation  of  the  efficiency, 
of  a  thermionic  diode  is 

J  V 

^  “  W+H  -H  +KaT  -  Q+X  (9) 

C  6l  J 

where  J  is  the  net  current  density,  is  the  output  voltage  (internal 
voltage  minus  the  potential  drops  in  leads  and  plasma),  W  is  the  net  thermal 
radiation  to  the  anode,  H  -H  is  the  net  energy  carried  by  the  electrons 

C 

from  the  cathode  to  the  anode,  KAT  is  the  heat  lost  by  conduction  through 
leads  and  plasma,  is  that  part  of  the  Joule  heating  in  leeuis  and  plasma 
that  returns  to  the  cathode,  and  X  is  the  support  losses.  Minor  terms,  such 
as  thermoelectric  effects,  have  been  neglected.  Space  charge  is  assumed  to 
be  neutralized,  either  intrinsicly  by  cesium  ionized  by  the  cathode  or 
by  some  other  means. 

This  expression  can  be  optimized  with  respect  to  both  lead 
resistance  and  load  resistance.  If  this  is  done,  the  values  of  efficiency 
are  functions  of  the  materials  of  construction  and  only  4  parameters  -  the 
cathode  temperature,  the  anode  temperature,  the  anode -to> cathode  spacing, 
and  the  cesium  bath  temperature.  A  complete  description  of  the  optimization 
procedure  is  presented  elsewhere. In  this  calculation,  several  effects 
ignored  by  most  other  authors  are  considered  and  shown  to  be  significant. 

In  particular,  the  effect  of  back  current,  thermal  conduction  through  the 
plasma,  and  electrical  resistivity  of  the  plasma  have  been  calculated. 

The  effects  of  varying  the  load  resistance  upon  the  behavior  of  both  the 
forward  and  back  currents  have  also  been  included. 
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The  diode  analyzed  for  this  paper  consists  of  two  plane  parallel 
tungsten  electrodes  separated  by  a  gap  filled  with  cesium  vapor.  The 
pressure  of  the  cesium  vapor  is  determined  by  the  temperature  of  the  cesium 
reservoirs.  The  temperatures  of  the  cathode  and  anode  were  considered  fixed 
by  means  of  temperature  reservoirs  and  are  taken  as  independent  parameters 
of  the  device.  The  thermal  emissivities  of  the  cathode  and  anode  were 
assumed  to  be  unaffected  by  cesium  coverage  and  were  taken  to  be  functions 
of  the  temperature  of  the  respective  surfaces  as  given  in  the  literature.'  ' 

The  work  functions  of  the  cathode  and  anode  were  determined  by  the  work  of 
Taylor  and  Langmuir  as  functions  of  the  appropriate  surface  temperature  and 
the  plasma  pressure. For  this  computation,  it  was  necessetry  to  fit  the 
data  of  Taylor  and  Langmuir  empirically  and  extrapolate  to  the  values  required. 
In  Figure  5  are  shown  curves  of  work  function  versus  surface  temperature  for 
several  cesium  bath  temperatures  as  calculated  from  this  empirical  formula. 
Iliese  curves  agreed  very  well  with  the  work  of  Taylor  and  Langmuir  within 
the  range  of  their  measurements.  The  values  of  plasma  thermal  conductivity 
and  plasma  resistance  as  functions  of  cesium  bath  temperature  were  taken 
from  the  experiments  of  Gottlieb  and  Zollweg.^^®^ 

These  calculations  were  carried  out  on  a  hi^  speed  digital  computer 
for  more  than  50OO  different  combinations  of  spacing,  and  cathode,  anode,  and 
cesium  bath  temperatures.  At  the  values  of  lead  and  load  resistance  which 
yielded  the  maximum  efficiency  for  each  of  the  given  combination  of  parameters, 
some  forty  different  properties,  such  as  power  delivered  to  load,  externeil 
voltage,  net  current,  etc.,  were  tabxilated  in  addition  to  the  efficiency. 

Some  sample  curves  of  efficiency  ere  given  in  Figs.  6-9.  These 
were  drawn  by  holding  two  of  the  four  parameters  constant  and  plotting  curves 
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of  constant  efficiency  as  a  function  of  the  remaining  two.  These  curves 
represent  the  efficiency  of  a  diode  which  has  been  optimized  with  respect 
to  both  lead  and  load  resistances.  In  Figure  6,  the  cathode  temperature 
and  spacing  have  been  fixed  and  the  efficiency  has  a  maximum  occurring  at 
a  unique  combination  of  cesium  and  anode  temperatures.  This  maximum  is  a 
general  occurrence  for  all  of  the  possible  combinations  of  spacing  and 
cathode  temperature,  although  its  precj.se  location  will  depend  upon  the 
particular  values  of  spacing  and  cathode  temperature  considered.  It  should 
be  mentioned  that  these  maxima  may  only  be  local  and  that  higher  efficiencies 
may  be  obtainable  at  lower  anode  temperatures.  However,  anode  temperatures 

less  than  600‘’K  were  not  considered  to  be  reasonable  from  a  design  standpoint. 

(3) 

A  complete  set  of  such  curves  have  been  presen  ed  elsewhere.  '  Figures  7  and 
8  are  drawn  similarly  to  Figure  6  except  the  temperature  fixed  was  that  of 
the  cesium  bath  or  anode,  respectively.  Curves  of  the  latter  two  types  do 
not  display  a  maximum  of  efficiency  since  higher  efficiencies  can  always  be 
obtained  by  going  to  a  sufficiently  higher  cathode  temperature. 

In  Figure  9>  'the  anode  and  cesium  temperatures  have  been  held 
constant  and  curves  of  constant  efficiency  have  been  drawn  as  a  function 
of  the  spacing  and  cathode  temperature.  For  a  fixed  value  of  cathode 
temperature,  it  can  be  shown  that  an  optimum  value  of  spacing  exists. 

However,  these  spacings  are  extremely  small  (less  than  2.5  microns)  and  not 
realizable.  Therefore,  for  spacings  larger  than  this,  it  can  be  seen  that 
efficiency  decreases  with  increasing  spacing. 

A  summary  of  the  maxima  obtainable  from  plots  of  the  type  of 
Figure  6  is  shown  in  Figures  10  and  11.  In  these  figures  are  shown  three 
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sets  of  curves  as  functions  of  cathode  temperature  -  maximum  attainable 
efficiency  and  the  anode  and  cesium  bath  temperatures  necessary  to  attain 
that  efficiency.  Each  set  consists  of  five  curves  corresponding  to  five 
different  spacings.  There  are  several  important  features  of  these  curves 
which  should  be  noted.  There  is,  for  a  fixed  spacing,  a  wide  range  of 
cathode  temperatures  in  which  the  ultimate  efficiency  is  relatively  indepen¬ 
dent  of  cathode  temperature .  Thus,  a  small  improvement  in  the  temperature 
capability  of  most  existing  materials  does  not  produce  a  large  improvement 
in  efficiency  of  operation.  Further,  the  anode  temperature  required  for 
maximum  efficiency  is  seen  to  be  within  the  range  900“K  to  1100®K  regardless 
of  the  cathode  temperature  or  spacing. 

It  should  be  emphasized  at  this  point  that  the  numerical  results 
presented  in  the  figures  and  in  the  above  discussion  refer  to  a  device  with 
tungsten  anode  and  cathode.  For  other  materisJ-s,  the  numerical  values  might 
change  considerably,  but  it  is  believed  that  the  qualitative  nature  of  the 
results  should  be  the  seune. 

It  must  also  be  recalled  that  these  calculations  are  based  on  the 
assumption  that  space  charge  has  been  eliminated.  Since  there  is  a  minimum 
cesium  temperature  at  which  this  is  true,  certain  areas  of  the  curves  presented 
(particularly  at  the  lower  cesium  temperatures)  are  Inaccurate.  Fortunately, 
from  the  limited  data  available  in  the  literature,  the  optimum  operating 
points  do  not  appear  to  be  in  this  excluded  area. 

IV.  Relationship  of  Efficiency  and  Lifetime 

With  the  aid  of  the  calculations  and  experiments  in  Parts  II  and 


III,  it  is  possible  to  determine  efficiency  and  lifetime  as  functions  of 
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only  four  parameters  —  cathode,  anode,  and  cesium  bath  temperatures  and 
the  cathode-anode  spacing.  If  this  determination  be  made,  then  for  a  fixed 
spacing,  a  maximum  efficiency  cein  be  found  for  a  given  minimum  evaporative 
lifetime.  The  measurement  of  mean  free  path  summarized  In  a  previous 
section  were  made  using  a  molybdenum  cathode.  The  relative  lifetimes 
deduced  from  this  value  should  not  be  seriously  different  for  tungsten, 
since  both  atomic  diameter  and  chemlcEd.  behavior  are  very  similar.  A  curve 
of  maxlotum  efficiency  as  a  function  of  the  evaporative  lifetime  for  tungsten 
Is  shown  In  Figure  12  for  a  123  micron  spacing  diode.  Similar  curves  can 
be  drawn  for  other  spaclngs.  Tills  efficiency  has  been  simultaneously 
maximized  with  respect  to  anode,  cathode,  and  cesium  temperatures  for  a 
fixed  lifetime.  In  principle,  this  procedure  could  be  continued  to 
determine  an  optimum  value  of  spacing.  However,  lanpractlcAbly  small  values 
of  spacing  result. 

The  extremely  large  values  of  lifetime  for  tungsten  shown  In  Fig. 
12  are  clearly  of  no  practical  significance.  However,  for  materials  of 
higher  vapor  pressure  than  tungsten,  an  analogous  curve  would  be  useful. 

In  particular,  ^  the  electron  emission  properties  of  molybdenum  In  the 
presence  of  cesium  were  the  same  as  those  of  tungsten,  then  the  lifetime 
scale  would  be  shifted  as  Indicated  In  Fig.  12.  It  Is  clear  that  under 
these  conditions,  a  required  lifetime  of  1  year  or  greater  will  result  In 
an  efficiency  less  than  20jt  for  molybdenum.  It  must  be  recsQ-led  that  this 
statement  Is  based  on  the  allowable  surface  recession  being  2.5  microns. 
However,  the  efficiency  In  the  region  of  Interest  for  molybdenum  Is  not  a 
sensitive  function  of  lifetime.  Allowing  10  times  the  recession  used  here 
does  not  markedly  increase  the  possible  efficiency. 
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If  accurate  emission  data  were  available,  calculations  of  the 
type  discussed  here  could  be  carried  out  for  other  materleds.  Similarly, 
lifetime  calculations  could  be  carried  out  If  measurements  of  mean  free  path 
for  other  materials  were  available.  While  experimented  Investigations  of  these 
data  cure  pres^'ntly  being  made  by  various  laboratories,  little  information  is 
available  In  che  literature. 

V.  Power  to  Weight  Ratio 

As  mentioned  In  Section  I,  it  may  often  be  more  desirable  to 

design  a  system  with  a  maxlmvuB  value  of  the  ratio  of  output  power  to  system 

weight  In  Ijeu  of  maximum  efficiency.  This  Is  an  extremely  difficult  problem 

to  solve  In  genered  since  the  weight  of  the  diode  Itself  Is  negligible,  and 

the  weight  of  the  system  Is  the  dominant  factor.  Ilhus,  the  construction  of 
the  system  must  be  known  In  great  detail.  For  exatople.  In  the  case  of  a 
reactor-powered  thermionic  converter,  a  reactor  core,  shield,  assembly, 
reflector,  and  radiator  would  have  to  be  designed  for  each  choice  of  spacing 
and  cathode,  anode,  and  cesium  tei^peratures  In  order  to  determine  the  weight 
of  the  system. 

It  Is  possible,  however,  to  say  a  few  things  In  a  semi -quantitative 
way  about  the  problem.  Consider  the  radiator,  which  Is  responsible  for  a  large 
part  of  the  weight  of  a  space  power  system.  Let  be  the  power  output,  P^ 
the  Input  power,  the  efficiency  of  conversion,  and  T  the  anode  temperature. 
Then  by  definition  of 


P  = 
o 


fio);: 


The  amount  of  power,  P^,  which  the  radiator  must  dump  as  waste  heats  is 
given  by 


^w  =  ^  -  ^o  = 

If  the  waste  heat  is  radiated 


P  =  C,AT 
w  la 


(12); 


where  A  is  the  surface  area  of  the  radiator  and  is  some  numeriCELl 
constant.  Further,  if  the  weight  of  the  radiator  is  proportional  to 
its  area. 


=  CgA  (15) 

When  combined,  the  above  equations  yield  for  the  ratio  of  output 
power  to  radiator  weight 


The  quantity  Ti(l-q)~^  is  a  monotonically  increasing  function  of 
T)  so  that  for  fixed  anode  temperature,  maximizing  the  efficiency  will 
maximize  as  well.  However,  in  many  cases,  the  anode  temperatures 

cannot  be  considered  fixed.  A  procedure  similar  to  that  used  in  the  preceding 
section  can  be  followed  using  the  maximum  power-to-weight  ratio,  instead  of 
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efficiency  as  the  criterion.  The  results  are  very  similar  to  those  of 
Section  III  in  that  higher  cathode  temperatures  result  in  higher  power-to- 
weight  ratios.  Figure  13  shows  the  results  of  optimizing  P for  a  5  mil 
spacing  tungsten  diode. 

A  comparison  of  the  resxilts  using  the  two  different  criteria, 
efficiency  (Fig.  10)  and  power- to-weight  ratio  (Fig.  13) >  shows  that  a 
significant  decrease  in  efficiency  is  the  price  that  must  be  paid  for  the 
maximum  power-to-weight  ratio.  This  is  particularly  true  in  the  region 
between  l600  and  2200“ K  cathode  teinperatures  where  the  optimum  efficiency 
(Fig.  10)  does  not  vary  rapidly  with  temperature. 

In  Figure  13  the  relative  values  of  the  power- to-wei^t  ratio  and 
the  optimum  anode  temperatures  are  also  plotted.  Ihe  peak  in  anode  tempera¬ 
tures  near  a  cathode  temperature  of  2500*K  is  almost  certainly  a  IocelL 
maximum,  since  the  power-to-weight  ratio  can  obviously  be  increased  by  going 
to  infinite  anode  temperatures,  using  a  "slightly  larger"  value  of  cathode 
temperature . 

VI.  Conclusions 

When  for  a  fixed  spacing  and  cathode  temperature  all  other  parameters 
are  optimized  using  efficiency  as  the  criterion,  the  resulting  optimum  efficiency 
is  a  relatively  insensitive  function  of  the  cathode  temperature  in  the 
temperature  range  l600  to  2200® K.  Clearly,  a  small  improvement  in  temperature 
capability  will  not  result  in  a  significant  improvement  in  efficiency.  Only 
a  large  Increase  in  operating  temperature  provides  a  marked  improvement.  In 
order  to  operate  a  thermionic  converter  of  the  type  discussed  here  at  25^ 
efficiency,  a  cathode  temperature  above  2200®K  is  necessary.  Only  tungsten 
is  capable  of  a  reasonable  lifetime  at  this  temperature. 


-38 


A  thermionic  convejter  using  tungsten  electrodes  operating  with 
a  cathode  temperature  of  l60C'  to  2200*K  would  have  an  efficiency  between 
15  and  20^  and  an  Indefinite  lifetime.  If  the  emission  characteristics  of 
molybdenum  are  similar  to  those  of  tungsten  a  molybdenum  diode  could  be 
constructed  with  a  lifetime  greater  than  one  year  at  an  efficiency  less 
than  2C$. 

The  introduction  of  cesium  into  a  thermionic  energy  converter 
has  two  opposing  effects  upon  the  efficiency.  Thermal  and  electrical 
losses  occur  in  the  cesium,  while  the  reduction  of  cathode  and  anode  work 
functions  allows  more  efficient  operation.  As  the  pressiire  of  cesium  is 
Increased  from  zero,  the  efficiency  is  increased  reaching  a  maximum.  Above 
this  critical  value  of  cesium  pressure,  the  efficiency  decreases  because 
of  the  increasing  electrical  and  thermal  losses  in  the  plasma.  The  attain¬ 
ment  of  maximum  efficiency  depends  critically  upon  establishing  the  optimum 
cesium  pressure.  The  use  of  a  cathode  material  possessing  a  larger. affinity 
for  cesium  would  make  possible  lower  cesium  pressure  to  obtain  a  given  work 
function.  Uie  plasma  losses  would  then  be  lower.  The  net  result  would  be 
a  higher  efficiency  of  operation.  In  short,  the  use  of  a  cathode  material 
of  higher  affinity  for  cesium  will  result  in  increased  efficiency.  It  can 
be  shown  that  the  converse  is  eilso  true.  Thus,  consideration  of  the  cesium 
affinity  is  of  considerable  importance  in  the  selection  of  a  cathode  material. 

From  examination  of  the  operating  characteristics  of  thermionic 
converters  (such  as  Figure  6)  it  is  seen  that  efficiency  is  a  relatively 
lnsensj.tlve  function  of  anode  and  cathode  temperature  and  of  spacing.  The 
cesium  pressure,  on  the  other  hand,  is  quite  a  critical  peurameter  and  should 
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be  controlled  very  closely.  The  lead  and  load  resistances  eire  also  critical 
(3) 

parameters.  Operation  away  from  these  crlticed.  values  can  decrease  the 
efficiency  by  a  sizeable  fraction.  Ihese  relationships  should  be  carefully 
considered  in  the  design  of  thermionic  converters. 

The  evaporative  lifetime  of  a  thermionic  converter  csjh  be  Increased 
by  increasing  the  spacing  and  the  cesium  pressure.  While  the  data  presented 
here  are  only  preliminary,  the  lifetime  of  a  cesium  filled  diode  will  probably 
be  at  least  10  times  that  estimated  by  the  use  of  evaporation  rates  into 
vacuum. 

The  data  necessary  for  calculations  of  the  type  described  are 
extremely  scarce  and  relatively  unreliable.  TOie  obtaining  of  more  extensive 
and  more  accurate  Information  would  allow  the  design  of  thermionic  converters 
to  be  made  in  a  relatively  routine  fashion.  For  tungsten  such  a  design  has 
been  made  and  the  device  is  presently  being  built. 

The  calcvilatlons  of  efficiency  have  been  checked  against  the  results 
of  the  few  experimental  devices  available  in  the  literature.  Where  sufficient 
information  is  given  to  edJLow  compsurison,  the  agreement  has  been  qiulte 
reasonable.  While  the  ceQ.culations  eure  not  to  be  considered  as  the  exact 
theoretlced.  results,  the  values  are  probably  quite  close  to  re8d.ity.  While 
the  model  of  plasma  is  not  entirely  realistic  and  the  effect  of  space  charge 
has  been  largely  Ignored,  the  errors  introduced  are  snail  in  the  regions  of 
interest.  Devices  of  35^  efficiency  are  clearly  Impracticail,  while  a  20)1 
efficient  device  can  probably  be  built  with  present  techniques  and  materials. 

Acknowledgements 

Some  of  the  work  described  herein  has  been  performed  under  Air 
Force  Contract  AF33(6l6)-d292> 


References 


4o 


1.  N.  S.  Baser,  J.  Appl.  Rliys.,  p.  l63  (i960); 

J.  M.  Hbuston,  J.  Appl.  Phys.  50,  p.  48l  (1959); 

J.  H.  Ingold,  J.  Appl.  Pbys.  3^7  p*  769  (196I); 

A.  Schock,  J.  Appl.  Phys.  p.  I56U  (196I) 

G.  N.  Hfttsopolous  and  J.  Kaye,  Proc.  I.R.E.  46,  p.  1574  (1958). 

2.  U.  S.  Patent  No.  2,9d0,dl9  Issued  to  G.  R.  Peaster,  April  Id,  196I. 

3.  L.  S.  Richardson  and  A.  E.  Fein,  mteria  Scientific  Report  No.  2, 
on  AP33( 616) -8292. 

4.  J.  V.  Edwards,  H.  L.  Johnston,  and  P.  E.  Blackburn,  J.  Am.  Chem.  Soc. 

I4,  p.  1539  (1952). 

R.  Speiser,  P.  E.  Blackburn,  and  H.  L.  Johnston,  J.  Electrochem.  Soc. 

106,  P.  52  (1959). 

J.  V.  Edvards,  H.  L.  Johnston,  and  P.  E.  Blackburn,  J.  Aa.  Chem.  Soc. 

p.  172  (1951). 

0.  Kubaschevski  and  E.  L.  Evans,  "Metallurgical  Thermochemistry",  Pergamon 
Press,  1958,  p.  335. 

5.  S.  Dushman,  "Vacuum  Techni<].ue",  Wiley,  1949;  P>  17* 

6.  Ibid,  p.  74 

7<  C.  V.  Maynard,  Nuclear  Sci.  and  Eng.  6,  p.  174  (1959)* 

8.  W.  E.  Forsythe  and  E.  M.  Watson,  J.  Qpt.  Soc.  of  Aaer.,  2^  p.  Il4  (1934). 
9<  J.  B.  Taylor  and  I.  lan^auir,  Phys.  Rev.  44,  p.  423  (1955) • 

10.  Unpublished  Research,  H.  Gottlieb  and  R.  ZoUweg. 


Lifetime  (seconds) 


4l 


CURVE  919066 


Fig.  1-Time  to  lose  2.5  microns  by  evaporation  from  molybdenum,  columbium, 
tantalum,  and  tungsten  as  a  function  of  temperature. 
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Fig.  2“Schematic  diagram  of  cell  for  determination 
of  evaporative  transport. 
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Fig.3-The  relative  transport  rate 
as  a  function  of  spacing. 
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4~Relativ6  lifetime  of  a  cathode  as  a  function  of  cesium  bath 
temoerature  and  soacina.  t- 
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Fig.6-Efficiency  os  o  function  of  anode  and  cesium 

bath  temperatures. 
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Fig.7- Efficiency  as  a  function  of 
cathode  and  anode  temperatures. 
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Fig.8- Efficiency  as  a  function  of  cathode 
and  cesium  temperatures 
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Fig.9- Efficiency  as  a  function  of 
spacing  and  cathode  temperature. 
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Fig.  1 1 -Optimum  onode  and  catium  temptroturti  oi  a  function  of 

cothodf  temptroturt. 
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Fig.  12-The  relationship  between  evaporative  lifetime  and  efficiency 
for  tungsten  and  molybdenum 
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THE  CALCUIATION  OP  THE  EFFICIENCY  OF  ENERGY  CONVERSION  BY 
THERMIONIC  EMISSION 


A.  E.  Fein  and  L.  S.  Richardson 


Introduction 

As  early  as  1953^  "the  possibility  of  using  thermionic  emission 
as  a  means  of  converting  thermal  energy  directly  into  electrical  energy 
was  discussed. Since  that  time,  a  number  of  papers  have  appeared  in 

the  literature  describing  the  operation  of  devices  utilizing  this  means 

(2\ 

of  energy  conversion. '  '  These  papers  have  used  the  efficiency  of  the 
device  as  a  criterion  of  comparison  with  other  energy  converters.  While 
many  of  the  authors  have  described  the  effects  of  varying  some  of  the 
many  parameters  governing  the  operation  of  a  thermionic  converter,  none 
have  considered  all  of  the  effects  simultaneously.  Certain  important 
effects,  such  as  back  current  and  varying  load  resistance,  have  not  been 
analyzed  carefully,  if  at  all.  In  this  paper,  an  attempt  is  made  to 
maximize  efficiency  by  simultaneously  optimizing  cesium  pressure  (including 
the  effects  on  work  function),  anode  temperature,  load  resistance  and  lead 
resistance.  The  effects  of  varying  cathode  temperature  and  anode-to-cathode 
spacing  are  also  included.  It  is  the  pxirpose  of  this  paper  to  construct  a 
model  predicting  the  behavior  of  a  thermionic  converter  for  the  use  of 
designers  or  experimenters.  Results  of  numerical  cong)utation  of  efficiency 
will  also  be  presented. 
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While  many  properties  are  of  importance  in  the  operation  of  a 
thermionic  converter,  once  the  materials  for  construction  of  cathode  and 
anode  have  been  chosen,  only  a  fevr  independent  parameters  remain.  The 
geometrical  configviration,  the  electrode  temperatures,  the  means  of  space 
charge  neutralization,  and  the  resistances  of  the  lead  and  load  are  the 
only  variables  which  can  be  adjusted.  The  most  common  means  of  space 
charge  neutralization  (recently  patented  by  G.  R.  Feaster^^^)  is  the  use 
of  a  cesium  atmosphere  in  the  interelectrode  gap.  Most  of  the  experimental 
cells  and  almost  all  of  the  converter  designs  use  geometry  which  is  equi¬ 
valent  to  plane  parallel  electrodes.  Since  the  lead  and  load  resistances 
can  be  simultaneo\jsly  optimized  to  yield  maximum  efficiency,  the  four  in¬ 
dependent  variables  in  a  conventionally  designed  thermionic  converter  are 
the  cathode  and  anode  temperatxires,  the  cesium  pressure,  and  the  inter¬ 
electrode  spacing.  A  later  section  will  show  that  the  anode  tenperature 
and  cesium  pressure  can  also  be  sim\iltaneously  optimized.  Thus,  once  the 
basic  design  of  a  converter  has  been  made,  the  choice  of  materials,  spacing, 
and  cathode  temperature  will  fix  all  other  parameters,  if  maximum  efficiency 
is  sought. 

II.  Description  of  the  Operation  of  a  Thermionic  Diode 

In  this  section,  a  mathematical  model  describing  the  behavior 
of  a  thermionic  diode  will  be  set  up  and  used  to  evaluate  the  efficiency 
of  such  a  device.  For  the  sake  of  mathematical  convenience,  the  operation 
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will  be  expressed  in  terms  of  variables  which  are  not  all  indei®ndent. 

(For  example,  thermal  emissivity  is  a  function  of  surface  temperatxore. ) 

This  will  not  affect  the  derivations  described  below.  Before  nmerical 
calculations  can  be  performed,  however,  the  actual  interrelationships  among 
the  physical  parameters  must  be  taken  into  account.  A  description  of  the 
manner  in  which  this  will  be  done  will  be  deferred  until  the  analysis  of 
the  operation  of  the  thermionic  diode  has  been  presented. 

A.  Description  of  Model  Used 

Stripped  of  the  numerous  parts  necessary  in  its  construction, 
the  thermionic  diode  can  be  basically  represented  by  the  schematic  of 
Fig.  1.  A  cathode  surface  in  contact  with  a  heat  reservoir  at  an  absolute 
temperature  thermally  emits  electrons  across  a  gap  to  an  anode  in  contact 
with  a  heat  reservoir  at  a  temperattire  T^.  The  electrons  are  returned  from 
the  anode  to  the  cathode  through  an  external  load. 

In  practice,  certain  effects  take  place  which  tend  to  prevent 
the  flow  of  a  reasonable  current.  The  most  serious  of  these  is  the 
formation  of  a  space  charge  in  the  space  between  the  cathode  and  anode. 

This  space  charge  forms  a  potential  barrier  for  the  electron  leaving  the 
cathode  and  can  reduce  the  n\anber  of  electrons  reaching  the  anode  (and 
therefore  the  current  through  the  load)  by  many  orders  of  magnitude. 

Two  different  methods  are  presently  in  use  to  overcome  the  problem  of 
space  charge.  The  simplest  of  the  two,  in  theory,  is  the  method  of  close- 
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(k) 

spacing.  It  is  known'  '  that  the  reduction  of  current  flow  by  space 
charge  is  very  sensitive  to  the  distance  between  the  cathode  and  anode. 

For  separations  of  0.0002  in.  or  less,  the  detrimental  effects  of  space 
charge  can  be  made  almost  negligible. ^  The  second  common  method  of 
reducing  or  eliminating  space  charge  is  by  the  use  of  a  postively  ionized 
vapor  (usually  cesim)  in  the  gap  between  the  anode  and  cathode.  The  electro¬ 
static  potential  set  up  by  the  ions  will  then  coimteract  the  retarding 
potential  of  the  electronic  space  charge  and  permit  an  Increased  flow 
of  electron  current. 

The  introduction  of  the  plasma  results  in  a  loss  of  energy 
available  to  the  load  by  conducting  heat  from  the  cathode  to  the  anode 
and  by  scattering  electrons  in  the  region  between  the  cathode  and  the 
anode.  In  addition,  the  plasma  may  modify  the  values  of  the  chemical 
potential  of  both  the  anode  and  cathode. 

In  this  paper,  it  will  be  assumed  that  the  space  charge  in  the 
region  between  the  anode  and  cathode  has  been  eliminated  by  using  a  plasma. 

It  will  be  assumed  that  the  sole  effect  of  the  plasma  on  the  behavior  of 
the  diode  is  to  introduce  a  "plasma  electrical  resistance,"  R^,  and  a 
"plasma  thermal  conductance,"  K^.  The  former  will  contribute  to  an  internal 
voltage  drop  inside  the  diode,  and  the  latter  to  a  thermal  conduction  of 
energy  across  the  cathode-anode  gap.  Also,  only  the  case  of  plane  parallel 
plates  one  square  centimeter  in  area  will  be  considered. 

In  this  regard,  it  must  be  remembered  that  the  value  of  resistance  mxist 
correspond  to  electrodes  of  \aiit  area,  and  that  net  current  densities 
calc\ilated  will  be  numerically  equal  to  currents. 
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Let  and  4  be  the  work  functions  of  the  cathode  and  anode 
c  a 

respectively.  (When  similar  qisantities  are  defined  for  both  the  cathode 

and  anode,  a  subscript  "c"  willrefer  to  the  former  and  "a"  to  the  latter. ) 

The  potential  seen  by  an  electron  in  going  from  the  cathode  to  the  anode 

will  then  be  represented  by  one  of  three  possible  cases  depicted  by  Fig. 

2.  The  mode  of  operation  will  be  classed  as  uphill  if  V  >  ♦  ~  9  (Fig. 

c  & 

2a),  downhill  if  V  <  ®  -  4  (Fig.  2b),  and  flat  if  V  =  4  -  4  (Fig.  2c) 

The  saturation  emission  current  densities  will  be  denoted  by 
J  and  J  .  By  Richardson's  equation 

C6  €IS 

-^CS  =  AT^exp  (-4^/KT^)  (l) 

•^as  =  AT^^exp  (-4^/^^)  (2) 

-2  -2 

where  A  is  Richardson's  constant,  equal  to  120  amp  cm  deg  K  j  and,  if 

k  is  Boltzmann's  constant,  and  e  the  electronic  charge,  K  =  k/e  =  8.6l6 

-5  -1 

x  10  volts  deg  K  .  In  the  flat  mode  the  net  current  density,  J,  from 
the  cathode  to  the  anode  is  simply  the  net  saturation  current  density,  J 

6 


J  = 


cs 


(3) 


In  general,  this  is  not  true,  and  the  net  current  must  be  written 


I 
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j  =  Jc  - 


where  in  the  \jphlll  mode 


Jc=  Jcs 


-(V  - 


J  =  J 
a  as 


and  In  the  downhill  mode 


(5) 


J  =  J 
c  cs 


J  =  J 
a  as 


exp  -  -  V  - 


(6) 


The  heat  input  to  the  cathode  and  power  output  to  the  load  must 

he  known  in  order  to  calculate  the  efficiency  of  operation  of  a  thermionic 

generator.  The  power  output  is  simply  the  load  voltage  times  the  current, 

where  the  load  voltage  is  equal  to  theemf,  V,  generated  by  the  diode  minus 

the  voltage  drops  due  to  plasma  resistance  (R  )  and  lead  resistance  (R. } . 

P  l 

fout  -  J(v  -  J  .  Rj,))  (7) 


By  the  first  law  of  thermodynamics  the  total  heat  input  to  the 
cathode  is  equal  to  the  total  energy  leaving  the  cathode.  The  latter  can 
be  expressed  as  the  sum  of  several  terms: 
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1)  Heat  radiated  from  cathode  to  anode,  W. 

2)  Support  and  end  losses,  X. 

3)  Heat  absorbed  by  electrons  leaving  the  cathode  and  reaching 
the  anode,  H  . 

4)  Heat  conducted  from  cathode  to  anode  by  lead  wires  and 
plasma,  H^. 

5)  Heat  retxirned  to  the  cathode  by  electrons  arriving  from  the 
anode,  H^^. 

6)  Heat  returned  to  cathode  from  the  Joule  heating  of  lead 
vires  and  plasma,  H,.> 

u 

The  first  four  effects  carry  heat  away  from  the  cathode  while  the  last  two 
bring  heat  to  the  cathode.  The  total  heat  input,  Hj^^,  to  the  cathode  can 
be  evaluated. 


Minor  losses  due  to  thermoelectric  effects  and  the  energy  required  to  ionize 
the  plasma  will  be  ignored. 


H^  =  W  +  H„  +  H  -  H  -  H_  +  X 
in  c  K  a  J 


(8) 


The  individual  terms  are  readily  calculable.  The  thermal  radiation  term 
W  is  given  by  the  Stefan-Boltzmann  law 


It  is  here  assumed  that  both  the  anode  and  the  cathode  are  gray  bodies. 
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W  =  OF  (9) 

where,  if  and  are  the  total  emisBivlties  of  a  parallel  cathode 
and  anode,  the  view  function  F 

F  =  (e^"^  +  -  1)“-^  (10) 

••12  •2 

and  0  is  the  Stefan-Boltzmann  constant  (5.67  x  10  watts  cm  deg  K  ). 
The  heat  carried  by  thermal  conduction  away  from  the  cathode  is 

where  and  are  the  thermal  conductances  of  the  plasma  and  the  leads 

respectively.  It  will  be  ass\Med  that  the  Joule  heat  returned  to  the 

cathode  is  simply  one-half  the  total  Joule  heat  produced  in  the  leads 
*  _ 

and  in  the  plasma.  Thus 


This  resvilt  is  obtained  rigorously  in  reference  2d. 


Hj  .  (1/2)  1  Bp)  ^ 


(12) 


In  order  to  compute  H  and  H  ,  the  mode  of  operation  must  be  considered. 

C  o 


From  Fig.  2,  it  is  seen  that  in  the  downhill  mode  an  electron  leaving  the 
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cathode  miast  have  absorbed  an  amount  of  heat  at  least  equal  to 
Similarly,  an  electron  arriving  at  the  cathode  from  the  anode  must 
also  have  absorbed  at  least  that  amount  of  heat.  For  the  uphill  mode, 

this  minimum  amovint  of  absorbed  heat  is  V  +  ®  .  Define 

a 

4  =  Max  (4^,  4^  +  V)  (13) 

The  value  of  4  is  thus  the  mlnimvon  heat  absorbed  by  an  electron  traversing 
the  diode  in  any  mode.  In  addition  to  this  minimum  heat  4,  there  is  some 
kinetic  energy  which  the  electron  may  possess.  If  the  electrons  immediately 
outside  the  cathode  or  anode  have  a  Maxwellian  distribution  characteristic 
of  the  temperature  of  the  surface  which  they  have  left,  it  can  be  shown  that, 
on  the  average,  this  additional  kinetic  energy  will  be  2KT  where  T  is  the 
temperature  of  the  emitting  sxtrface.  Thus 

(4  +  2KT^)  (14) 

and 

=  Jg  (4  +  2KTg)  (15) 

The  efficiency,  defined  as  the  ratio  of  output  power  to  heat  in,  is  then 
given  by 
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T1  = 


J  [V  -  J  (Rj^  +  Rp)] 


W  +  J„(«+2KP„)  +  (Kj^+Kp)(T^-T^)  -  J,(<t+2KT,)  -  j2(R^+R^)/2  +  X 


(16) 


L-*p^ 


In  practice,  the  load  resistance,  1®  ®  quantity  which 

can  he  varied  at  will.  However,  mathematically  the  internal  voltage 


V  =  J(R^  .  Ej,  t 


(17) 


is  much  more  convenient  to  use.  Therefore,  in  lieu  of  R^,  Rj^  and  R2.oad 

being  considered  as  three  independent  variables,  R^,  R^^  and  V  will  be 

used,  the  transformation  being  non-singular.  It  must  be  ke^in  mind, 

therefore,  in  the  following  analysis  that,  i^J  for  example,  a  variation 

of  Rj^  at  constant  V  is  considered,  a  variation  in  R^  and/or  Is 

necessary.  The  value  of  the  variable  V  is  not  completely  optional.  There 

are  limits  on  V,  and  values  outside  of  these  limits  are  not  physically 

obtainable.  It  is  only  the  load  resistance  which  can  be  physically  varied 

from  zero  to  infinity,  short  and  open  circviit  respectively.  Each  case 

yields  an  equation  which  may  be  solved  for  the  internal  voltage.  The 

minimum  value  for  V,  V  ,  is  obtained  from  the  short  circuit  equation 
•  sc' 


0  .  V  -  J  (Ej,  .  Ep)  at  V  . 


(18) 


Equation  (l8)  is  transcendental  and  cannot  be  solved  analytically.  The 
open  circuit  voltage  is  determined  by  the  condition 
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J  =  0  at  V  =  V  (19) 

oc  ' 

An  analytic  solution  to  (19)  for  V  can  be  obtained  from  (4),  (5)  and 

OC 

(6) 


oc  c  a  '  cs'  as 


(20) 


where 


T  =  if  J  „ 
c  cs 


>  J  and  T  =  T.  if  J  .  >  J 


as 


as 


cs 


B.  Optimization  of  ^ad  Resistance 

The  efficiency  (16)  can  be  optimized  with  respect  to  the  lead 
resistance  at  fixed  V  as  follows: 

If  L,  p,  and K are  respectively  the  cross-sectional  area, 
length,  electrical  resistivity,  and  thermal  conductivity  of  the  leads 


\  -  PL/A,. 

(21) 

=  kAc/L 

(22) 

By  the  Wiedemann-Franz  law  p  and  k  of  a  metallic  conductor  at  a  temperature 
are  related  by  the  equation 
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P  =  X  =  2.44  X  10‘®  watt-ohm-deg  K"^  (23) 

If  T  is  asstmed  to  be  the  average  of  T  and  T 
0  c  a 

=  (\/2  R^)  (T^  +  Tg)  (24) 

Thus  can  be  replaced  in  (l6)  by  its  value  in  terms  of  R^^ 
from  (24)  and  obtain  an  expression  which  can  be  reduced  to 


n  =  Plo 


(25) 


where  is  the  efficiency  of  an  identical  diode  with  ideal  leads  (i.e., 

Rj^  =  0,  =  0) 


J[V  -  J^] 


W+J  (<I>+2KT  )  -  J  ('In-2KT„).+  K  (T  -T„)  -  /2  +  X 

c'  c'  a'  a'  p'  c  a'  p' 


(26) 


and 


P  -  (1-7 )  "  2 

G  =  (V-JRp)^  [(T^^-  T^2)]  ^ 


(27) 

(28) 


(29) 
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From  (17)  therefore 


y 


+  R 


load 


(50) 


Thiis  physically  the  parameter  y  is  the  fraction  of  the  external  emf  of 
the  diode  which  appears  across  the  leads)  or  equivalently,  it  is  the 
ratio  of  the  lead  resistance  to  the  total  external  resistance.  Its  range 
therefore  is  from  0  to  1. 

The  efficiency  can  then,  in  principle,  be  maximized  with  respect 
to  the  lead  resistance  Using  the  representation  of  (25)  and  notltiig 
that  R^  enters  only  through  the  dimensionless  parameter  7,  it  is  found 
that  the  efficiency  optimized  with  respect  to  R^^  is  (qiantities  evaluated 
at  that  value  of  R_  giving  maximum  efficiency  will  be  denoted  with  an 
overhead  bar) 


n 


(51) 


where 


5=^7^- " 

(52) 

v/1  +  B  +  1  -  q 
'0 

\  =  [(V  -  JHp)/J][('/l  +  B  -  1)/B] 

(55) 

7  =  Jgj^/(V  -  JRp). 

(54) 

68 


®  (55) 

As  the  value  of  B  ranges  from  zero  to  infinity,  the  value  of 
y  will  range  from  1/2  to  zero.  Thus  under  no  conditions  will  more  than 
half  the  external  emf  appear  across  the  leads  when  they  are  optimized. 

Obviously,  3  <  1.  The  quantity  p  is  the  fraction  to  which  the 
efficiency  of  the  ideal  diode  with  ideal  leads  (zero  electrical  resistance 
and  zero  thermal  conductivity)  has  been  reduced  when  the  physical  leads 
have  been  optimized.  Considered  as  a  function  of  alone  (v,  T^,  T^ 
fixed),  3  can  be  shown  to  be  monotonically  decreasing.  Thus  the  tendency 
of  the  leads  to  decrease  the  efficiency  is  greater  for  an  ideal  diode  with 
a  large  efficiency.  This  can  be  ecplained  qualitatively  from  the  fact  that 
the  higher  efficiency  diodes  have  in  general  the  larger  electric  currents 
and  that  the  larger  the  current,  the  larger  the  voltage  drop  in  the  leads. 
However,  rj  is  still  a  monotonically  increasing  function  of 

If  non-optimum  leads  are  used,  a  decrease  in  efficiency  will 

result.  This  effect  will  be  considered  for  the  most  pessimistic  case 

by  adjusting  the  parameters  so  that  B  is  as  small  as  possible.  This 

can  be  done  by  making  G  as  small  as  possible.  For  the  values  =  0°K, 

T  =  3000°K,  and  V  -  JR  =1  volt,  G  is  equal  to  4.55*  Fig.  3  shows  the 
c  p 

results  of  calculations  of  3  as  a  function  of  and  y.  Also  shown  in  that 
figure  is  the  locus  of  the  maximum  of  the  curves,  3.  Several  features  of 
these  curves  are  of  interest.  At  a  critical  value  of  y  equal  to 
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(0.4T  in  this  case),  the  curves  for  different  values  of  intersect.  For 
7  greater  than  this  critical  value,  p  is  a  monotonically  decreasing  function 
of  hut  very  weakly  dependent  \3Son  For  7  less  than  this  critical 
value,  p  is  a  strong  monotonically  Increasing  function  of 

Curves  of  y,  rj,  and  p  as  functions  of  are  presented  in  Fig. 

4  for  the  values  of  Tg,  and  V  used  above. 

It  is  clear  from  Figs.  3  and  4  that  there  is  a  serious  loss  in 
efficiency  due  to  the  leads.  Even  when  this  loss  is  minimized  by  opti¬ 
mizing  the  value  of  the  lead  resistance,  it  is  seen  that  for  the  case 
considered,  the  ideal  efficiency,  is  decreased  by  amounts  ranging 
from  1^  (at  t|^  =  0.1)  to  42)b  (at  tJq  =  1.0).  For  less  pessimistic  cases, 
the  loss  is  less  severe,  but  still  appreciable.  It  should  be  enqphasized 
that  the  ctirves  of  Fig.  are  for  optimized  leads.  For  non-optimized  leads, 
the  loss  of  efficiency  is  much  greater  as  can  he  seen  from  Fig.  5- 

It  is  seen  from  Fig.  4-  that,  for  almost  the  entire  range  t)  is, 
to  a  good  approximation,  a  linear  faction  of  This  fact  is  borne  out 

hy  expanding  equation  (31)  in  a  series  about  =  1. 


C(1  - 
0.1 


(36) 


For  the  case  considered,  (G  =  55)> 


1  -  0.82  (1  -  tJq)  -  O.Ol+T  (1  -  +  ••• 


ii  =  .57»+ 


(37) 
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The  high  degree  of  linearlity  of  t]  as  a  function  of  observed 
for  this  case  does  in  fact  occur  in  general.  This  is  seen  by  examining 
the  coefficients  of  the  linear  and  second  order  terms  in  (56).  For  the 
case  treated,  the  ratio  of  the  coefficient  of  the  second  order  term  to 
that  of  the  first  order  term  is  0.05?.  For  values  of  G  >  h. 55,  the  ratio 
decreases  vmtil  at  G  =  8,  it  becomes  zero.  For  larger  values  of  G,  the 
ratio  becomes  negative  and  obtains  a  maxinium  magnitude  of  0.042  at  G  *  35* 
Thus  the  special  case  which  has  been  considered  actually  exhibits  the 
largest  deviation  from  linearity.  Therefore,  to  a  good  approximation 
(less  than  %  error)  tj  (for  >  0.1)  can  be  computed  from  only  the 
first  two  terms  of  (36).  Similarly,  for  small  equation  (51)  can 
be  e:q)anded  in  a  series  of  powers  of  Thus  Tj  can  be  written 


{  (Vl  +  G  -  1)(1  +  G)"^/^ 


1  +  GtIo  + 


n  = 


1  -  2^^J2Q  +  2(t,^/2G)  ... 


T].  >  0.1 


n  <  0.1 


(58) 


For  the  lower  series  of  (58),  it  is  again  noted  that  for  the  case  of 
minimum  G  treated,  the  fit  is  the  poorest,  and  that  for  larger  values 
of  G  the  approximation  improves. 
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C.  Optimization  of  Load  Resistance  (internal  Voltage) 

The  assmption  has  often  been  made  in  the  literature  that  the 
flat  mode  of  operation  yields  the  maximum  attainable  efficiency  for  a 
given  diode.  This  assvmqptlon  is  not  necessarily  valid,  as  a  simple 
numerical  example  will  show.  If  R  =  o,  *  =  2-5  volts,  *  =  1,5  volts, 

=  2200°K,  Tg  =  0°K,  and  and  =  .25,  then  =  34. 5^t  when  operating 
in  the  flat  mode.  Using  the  same  parameters,  but  operating  in  the  uphill 
mode,  so  that  V  =  I.5  volts,  Is  equal  to  41.5^*  At  the  same  time,  a 

P  p 

reduction  of  power  density  from  1090  watts /cm  to  II7  watts /cm  and  a 

O  P 

reduction  of  current  density  from  1090  amps /cm  to  77.8  amps /cm  takes 
place.  Thvus  not  only  has  the  efficiency  increased  upon  increasing  the 
output  voltage,  but  both  the  current  density  and  the  power  density  have 
been  considerably  reduced.  The  reduction  of  the  power  density  is  of 
significant  value  because  a  comparable  power  must  be  dissipated  at  the 
anode,  and  the  dissipation  of  large  quantities  of  power  can  pose  a  serious 
engineering  problem.  Similarly,  the  handling  of  large  current  densities  is 
difficviltj  thvis  the  reduction  of  the  output  current  density  is  also  of 
engineering  importance. 

In  this  exanqjle,  the  decrease  of  the  current  and  power  densities 
was  accompanied  by  a  significant  Increase  in  efficiency.  Cases  can  be 
constructed  in  which  the  change  in  efficiency  is  not  significant,  but 
in  which  operation  in  the  uphill  mode  decreases  the  current  and  power 
densities  by  appreciable  amounts. 


I 

I 
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The  reason  for  this  Increase  in  efficiency  can  he  seen 
qualitatively  by  considering  the  simple  case  of  ideal  leads,  no  plasma 
or  support  losses,  and  no  back  CTirrent.  Then  (l6)  can  be  written 


As  long  as  the  thermal  radiation,  W,  is  small  compared  to  J^,  the  efficiency 
will  increase  with  increasing  V.  (*  =  V  +  it  in  the  uphill  mode).  Event\ially, 

ci 

J  becomes  so  small  that  W/J  is  no  Icnger  negligible.  At  this  point,  the 
c  c 

efficiency  begins  to  decrease  with  increasing  V.  When  a  maximum  efficiency 
occurs  in  the  uphill  mode,  it  can  be  shown  that  this  efficiency  is  independent 
of  the  cathode  work  function,  «  . 

Nmerlcal  calculations  of  (l6)  show  that  this  effect  exists  as 
well  for  that  more  accurate  case.  Further,  it  is  possible  for  the  optimum 
efficiency  to  occur  in  the  downhill  mode.  The  fact  that  this  can  occur  can 
easily  be  explained  by  considering  the  case  of  a  diode  whose  parameters  are 
such  that  J..  <  J  Thus  the  only  possible  mode  of  operation  is  the  down- 

CS  Q8 

hill  mode  and  the  maximm  efficiency  must  occur  in  that  mode. 

The  efficiency  of  (l6)  can  be  simultaneously  maximized  with 
respect  to  V  and  by  setting  equal  to  zero  the  partial  derivatives  of 
T)  with  respect  to  these  variables.  This  procedure  results  in  the  equation 


for  V 
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J  +  [V  -  2J  +  Rp)]  -  ar 
dj 

-  ^  +  2KT^)  W  (R^  +  Rp)]  =  0 


(**0) 


where  tj  and  are  given  as  functions  of  V  by  (Jl)  ff.  Equation  (40)  is 
actually  two  equations  —  one  for  the  downhill  and  one  for  the  uphill 
mode  of  operation.  It  can  be  shown  that  each  of  these  two  equations 
always  has  one  and  only  one  solution  for  V  between  V  and  V  .  If  the 
solution  to  the  equation  for  the  downhill  anode  lies  between  V  and 

8C 

there  is  a  maximum  for  the  efficiency  in  that  mode.  Similarly, 
if  the  equation  for  the  uphill  mode  has  a  solution  between  and 

there  is  a  maximum  in  the  uphill  mode.  If  neither  of  these  two  conditions 
result,  the  maximum  efficiency  results  in  the  flat  mode  (v  =  -  ♦g). 


In  principle,  it  is  possible  for  a  maxlmian  to  occur  in  each  mode. 
However,  the  range  of  parameters  for  which  this  can  happen  is  small. 

In  several  thousand  cases  calc\ilated,  this  situation  has  never  occ\irred. 


Once  the  value  of  V  resulting  in  the  maximum  efficiency  is  found, 
the  value  of  Rj^  necessary  for  it  to  be  obtained  is  given  by  (53)-  And 
from  these  two  quantities  all  other  properties  at  this  optimum  point  can 
be  found  (e.g.,  power  out,  net  current,  etc.). 
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The  model  of  the  diode  used  for  the  above  calculation 
possesses  several  limitations  with  respect  to  the  manner  in  which  the 
plasma  is  included.  These  are  important  enough  to  warrant  further  dis¬ 
cussion.  It  has  been  assmed  that  "sufficient"  plasma  has  been  intro¬ 
duced  so  that  the  space  charge  is  eliminated.  Although  ejjperlments  with 
plasma  filled  diodes  seem  to  indicate  this  can  be  done,  only  the  crudest 
estimates  are  available  of  the  conditions  under  which  this  is  possible. 
Further,  it  is  known  that  the  Introduction  of  plasma  can  result  in 
potential  sheaths  at  the  electrodes.  In  lieu  of  considering  this 
difficult  and  as  yet  unsolved  problem,  the  procediare  described  assumes 
that  the  plasma  can  be  considered  ohmic  and  that  an  effective  electrical 
resistance  can  be  found.  Limited  experimental  data  indicates  that  this 
assumption  is  not  too  inaccvirate.  Finally,  the  possibility  of  an  arc 
being  established  is  ignored.  Although  it  is  definitely  known  that  an 
arcing  mode  can  be  set  up,  this  possibility  does  not  invalidate  the 
above  calculations.  What  it  means,  however,  is  that  if  the  solution  for 
the  operating  point  found  by  the  described  procedure  falls  in  a  region  of 
operation  which  will  support  an  arc,  it  is  not  a  valid  solution  for  the 
maxintum  efficiency.  The  maximum  efficiency  in  the  region  where  no  arc 
occurs  miist  then  be  found.  In  either  case,  to  determine  the  optimm 
efficiency,  the  operation  in  the  arcing  mode  must  be  analyzed  and  the 
maximum  efficiency  obtained  there  con5)ared  to  the  maximum  efficiency  of 
normal  operation  already  described.  Again  a  treatment  of  the  operation 
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of  a  diode  in  the  presence  of  an  arc  is  a  difficult  unsolved  problem. 

By  comparing  the  optimum  operating  points  calculated  by  the  formulas 
given  here  to  the  limited  e3q>erimental  data  available  on  arcs^  it  is 
found  that,  in  general,  these  cosgputed  optimum  points  should  not  sustain 
arcs. 

D.  Empirical  Evaluation  of  Parameters 

In  order  to  be  able  to  calculate  the  quantities  dldcussed,  it 
is  necessary  to  know  the  values  of  many  parameters  such  as  work  functions 
and  emissivlties.  In  reality,  quantities  such  as  these  are  not  independent, 
but  depend  on  the  pressure  of  cesiiam  and  the  temperatures  of  the  electrodes. 
This  section  is  devoted  to  describing  the  relation  which  exists  eunong  these 
quantities . 

The  relationship  of  effective  work  function  to  cesium  pressure 

and  sxirface  temperature  has  not  been  well  determined.  The  pioneering  work 

(l) 

of  Langmuir' ' '  provided  the  basis  for  the  present  empirical  equation. 

While  no  theoretical  justification  of  the  equation  can  be  given,  the 
results  are  in  close  agreement  with  the  original  data  and  with  the  extra¬ 
polations  of  other  investigators. 

The  effective  work  fmction  of  cesium-coated  tungsten  was  assumed 
to  be  a  function  of  only  the  fractional  coverage,  6.  The  data  of  Langmuir 

(7^ 

and  Taylor'"  are  plotted  in  Fig.  5.  The  empirical  equation 
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<b  =  4.52  -  8.59  e  +  6.3  ©  <  .86 

t  =  1.81  ,  ©  >  .86 


(41) 


was  fitted  to  this  data.  The  resuilts  of  this  equation  are  also  plotted 
in  Fig.  5* 

The  estimation  of  the  fractional  coverage  as  a  function  of 
cesivun  pressiare  and  surface  temperature  was  also  made  by  the  vise  of  a 
set  of  empirical  equations.  Fig.  6  shows  the  relation  between  surface 

(7) 

temperatxrre  and  coverage  for  two  particular  cesium  temperatures.'*' 

This  can  be  well  fitted  by  the  equation 


©  =  ri  +(T/T^)'^] 


“1 


(42) 


where  T  and  g  are  constants,  both  of  which  are  functions  of  the  cesium 
o  ' 

pressure  (or  arrival  rate). 

If  the  relationship  of  cesium  bath  temperature^  Tq^,  to  cesivun 

/qA 

pressiire,  P,  given  by  Nottingham'  ^  is  used^ 


=  2.45  X  10"®  e 


(43) 


where  Tj,^  is  in  deg  K  and  P  in  ram  Hg.  The  relationship  of  c  and  T^  to 
the  tenrperature  of  the  cesium  bath  are  found  to  be 
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c  =  2.73  +  720/T^g 

1  .1; 

0.385  (10  +  1/Tcs) 

o 

Figure  6  shows  the  resiilts  of  using  these  empirical  equations  for  the 

two  temperatures  shown.  Figure  7  shows  the  relationship  of  work  function 

to  cesixim  temperature  and  siarface  tempera tvire .  This  is  in  fair  agreement 

with  the  res\ilts  of  Schock'  especially  if  the  cesium  temperature  listed 

by  Schock  rather  than  the  pressure  is  used  to  compare  values.  One  particular 

difference  js  of  interest.  The  meastarements  of  Langmuir  indicate  that  a 

minimum  work  function  is  obtained  when  9  is  equal  to  2/5.  When  a  full 

monolayer  is  present,  the  work  function  is  assumed  to  be  that  of  cesium. 

This  higher  value  of  the  work  fvinction  at  full  coverage  has  a  serious 

effect' on  the  efficiency  when  low  anode  temperatvures  are  used. 

The  thermal  emissivities  of  tungsten  is  given  by  several 

references  as  a  function  of  xemperatiire.  The  data  used  here  are  those 

(9) 

of  Forsythe  and  Watson. 

Values  of  plasma  electrical  resistance  and  thermal  conductance 
are  not  available  in  the  literature.  Gottlieb  and  ZoUweg^^^^  have 
recently  meas\ired  the  thermal  conductance  of  cesium  plasma.  An  empirical 
fit  to  their  results  is  shown  in  Figure  8.  These  data  were  fitted  well 


by  the  empirical  equation. 
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Kp  =  1.5  X  10"^  |l  -  exp  (-1+0  Pd)|  (1+5) 

where  P  is  the  cesitan  presstire  in  mm  of  Hg  (by  (43)),  d  the  spacing 
in  cm  and  Kp  the  thermal  conductivity  of  the  plasma  in  watt-cm/deg  K. 

The  electrical  resistance  was  estimated  from  J-V  curves  taken  by 
Gottlieb  and  Zollweg.  The  resvLLts  are  quite  preliminary  in  natvire. 

For  the  calculations  discussed  hereafter,  the  plasma  resistivity,  pp, 
was  estimated  by  the  equation 

Pp  =  2.67P  (46) 

where  P  is  the  cesim  pressure  in  mm  Hg.  (by  (45))  and  pp  is  in  ohm  cm. 

III.  ResTolts  of  Numerical  Calculations 

Because  empirical  data  concerning  the  effect  of  cesium  coverage 
on  the  thermionic  properties  of  materials  is  well  known  only  for  tungsten, 
calcvilations  were  made  only  for  a  device  in  which  both  the  anode  and  cathode 
were  made  of  tungsten. 

Calculations  of  the  optinrum  operating  point  according  to  the 
method  described  above  were  carried  out  on  a  high  speed  digital  ccMputer 
for  more  than  5000  different  combinations  of  spacing,  and  cathode,  anode, 
and  cesium  bath  temperatures.  At  the  values  of  lead  and  load  resistance 
which  yielded  the  maximum  efficiency  for  each  of  the  given  combination  of 
parameters,  some  forty  different  properties,  such  as  power  delivered  to 
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load^  external  voltage,  net  current,  etc.,  were  tabulated  in  addition  to 
the  efficiency  (a  sainple  of  such  a  table  is  shown  in  Table  I).  As  already 
shown,  for  given  materials  of  construction,  the  efficiency  t),  optimized  with 
respect  to  load  and  lead  resistance  is  a  function  of  only  four  variables  - 
cathode-anode  spacing,  and  the  tei^perature  of  the  cathode,  anode,  and 
cesivim  reservoir.  Because  of  the  nmber  of  independent  variables,  it  is 
possible  to  plot  lines  of  constant  efficiency  in  a  large  number  of  ways. 

By  keeping  two  of  the  parameters  constant,  efficiency  contours  can  be 
plotted  as  a  function  of  the  other  two.  Figures  9,  10,  and  11  show 
respectively  sample  results  of  keeping  fixed  cesim  temperatiire  and 
spacing,  anode  temperature  and  spacing,  and  cesium  temperature  and  anode 
temperature.  It  was  found,  however,  that  the  most  useful  type  of  plot  was 
obtained  by  fixing  the  spacing  and  cathode  temperature  and  drawing  curves 
of  constant  efficiency  as  a  function  of  anode  and  cesiian  reservoir  tempera¬ 
tures.  A  set  of  60  (5  spacings  x  12  cathode  temperatures)  of  such  plots 
are  shown  in  Figure  12.  It  should  be  noticed  that  in  each  case,  the 
efficiency  has  a  maximum  occxirring  at  amique  combination  of  cesium  and 
anode  temperatures.  This  maximum  is  a  general  occurrence  for  all  of  the 
possible  combinations  of  spacing  and  cathode  temperature,  althou^its 
precise  location  will  depend  upon  the  particular  values  of  spacing  and 
cathode  temperature  considered.  It  shoiild  be  mentioned  that  these  maxima 
may  only  be  local  and  that  higher  efficiencies  may  be  obtainable  at  lower 
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anode  tenrperatures .  However,  anode  temperatures  less  than  600°K  were  not 
considered  to  be  reasonable  from  a  design  standpoint.  Further,  the  anode 
teinperatvire  can  never  be  less  than  the  temperature  of  the  cesium  reservoir, 
for  the  latter  must  be  the  lowest  temperature  in  the  diode  so  that  it  will 
determine  the  pressure  of  the  cesium  vapor.  Curves  of  the  type  in  Figs.  9, 

10,  and  11  do  not  display  a  maximum  of  efficiency  since  hlgho'  efficiencies  can 
always  be  obtained  by  going  to  a  stiff Iclently  higher  cathode  temperatxire. 

If  Fig.  11  is  considered  at  a  fixed  value  of  cathode  temperature,  it  can 
be  shown  that  an  optimum  value  of  spacing  exists.  However,  these  spacings 
are  extremely  small  (less  than  2.5  microns)  and  not  realizable.  Therefore, 
for  spacings  larger  than  this,  the  efficiency  decreases  with  increasing 
spacing. 

A  summary  of  the  maxima  obtainable  from  plots  of  the  type  of 
Fig.  12  is  shown  in  Figs.  15  and  14.  In  these  figures  are  shown  three 
sets  of  curves  as  functions  of  cathode  temperature  -  maximum  attainable 
efficiency  and  the  anode  and  cesium  bath  tenperatures  necessary  to  attain 
that  efficiency.  Each  set  consists  of  five  curves  corresponding  to  five 
different  spacings.  There  are  several  important  features  of  these  cvurves 
which  shouldbe  noted.  There  is,  for  a  fixed  spacing,  a  wide  range  of 
cathode  temperatures  in  which  the  ultimate  efficiency  is  relatively  in¬ 
dependent  of  cathode  tenperatures.  Thus,  a  small  improvement  in  the 
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ten®erat\ire  capability  of  most  existing  materials  does  not  produce  a  large 
improvement  in  efficiency  of  operation.  Further,  the  anode  temperature 
required  for  maximum  efficiency  is  seen  to  be  within  the  range  900°K 
to  1100^  regardless  of  the  cathode  temperature  or  spacing. 

It  should  be  reemphasized  at  this  point  that  the  numerical  results 
presented  in  the  figures  and  in  the  above  discussion  refer  to  a  device  with 
tungsten  anode  and  cathode.  For  other  materials,  the  numerical  values  ml£^t 
change  considerably,  but  it  is  believed  that  the  qualitative  nature  of  the 
results  should  be  the  sama 

It  must  also  be  recalled  that  these  calculations  are  based  on 
the  assumption  that  space  charge  has  been  eliminated.  Since  there  is 
a  minimum  ceslvnn  temperature  at  which  this  is  true,  certain  areas  of  the 
curves  presented  (particularly  at  the  lower  cesium  temperatures)  are  in¬ 
accurate.  Fortunately,  from  the  limited  data  available  in  the  literature, 
the  optimum  operating  points  do  not  appear  to  be  in  this  excluded  area. 

If  accurate  emission  data  were  available,  calculations  of  the 
type  discussed  here  could  be  carried  out  for  other  materials.  While 
experimental  investigations  of  these  data  are  presently  being  made  by 
various  laboratories,  little  Information  is  available  in  the  literature. 

IV.  Conclusions  and  Recommendations  for  Further  Study 

The  following  are  scmie  conclxisions  we  have  drawn  from  this 
study  and  some  recommendations  for  further  work. 
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When  for  a  fixed  spacing  and  cathode  temperature  all  other 
paremeters  are  optimized  vising  efficiency  as  thf  criterion,  the  resulting 
optimum  efficiency  is  a  relatively  insensitive  function  of  the  cathode 
temperature  in  the  temperature  range  l600  to  2200°K.  Clearly,  a  small 
improvement  in  temperatvire  capability  will  not  result  in  a  significant 
improvement  in  efficiency.  Only  a  large  increase  in  operating  temperature 
provides  a  marked  Improvement.  In  order  to  operate  a  thermionic  converter 
of  the  typ;?  discussed  here  at  2%  efficiency,  a  cathode  temperature  above 
2200°K  is  necessary.  A  thermionic  converter  using  tungsten  electrodes 
operating  with  a  cathode  temperature  of  1600  to  2200°K  would  have  an 
efficiency  between  15  and  20^. 

The  introduction  of  cesivan  into  a  thermionic  energy  converter 
has  two  opposing  effects  upon  the  efficiency.  Thermal  and  electrical 
losses  occur  in  the  cesivnn,  while  the  reduction  of  cathode  and  anode 
work  functions  allows  more  efficient  operation.  As  the  pressure  of 
cesium  is  increased  from  zero,  the  efficiency  is  Increased  reaching  a 
maximvim.  Above  this  critical  value  of  cesium  pressvire,  the  efficiency 
decreases  because  of  the  increasing  electrical  and  thermal  losses  in  the 
plasma.  The  attainment  of  maximvim  efficiency  depends  critically  upon 
establishing  the  optimum  cesium  pressure.  The  use  of  a  cathode  material 
possessing  a  larger  affinity  for  cesivira  would  make  possible  lower  cesivan 
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pressure  to  obtain  a  given  work  function.  The  plasma  losses  would  then 
be  lower.  The  net  result  would  be  a  higher  efficiency  of  operation.  In 
short,  the  use  of  a  cathode  material  of  higher  affinity  for  cesium  will 
resvilt  in  increased  efficiency.  It  can  be  shown  that  the  converse  is 
also  true.  Thus,  consideration  of  the  cesium  affinity  is  of  considerable 
Importance  in  the  selection  of  a  cathode  material. 

From  examination  of  the  operating  characteristics  of  thermionic 
converters  (such  as  Fig.  12),  it  is  seen  that  efficiency  is  a  relatively 
insensitive  function  of  anode  and  cathode  temperature  and  of  spacing. 

The  cesltm  pressure,  on  the  other  hand,  is  quite  a  critical  parameter  and 
sho\ild  be  controlled  very  closely.  The  lead  and  load  resistances  are  also 
critical  parameters.  Operation  away  from  these  critical  values  can 
decrease  the  efficiency  by  a  sizeable  fraction.  These  relationships  should 
be  carefully  considered  in  the  design  of  thermionic  converters. 

The  calcvilations  of  efficiency  have  been  checked  against  the 
results  of  the  few  experimental  devices  available  in  the  literature.  Where 
sufficient  information  is  given  to  allow  conparison,  the  agreement  has  been 
quite  reasonable.  While  the  calculations  are  not  to  be  considered  as  exact 
theoretical  results,  the  values  are  probably  quite  close  to  reality.  While 
the  model  of  the  plasma  is  not  entirely  realistic  and  the  effect  of  space 
charge  has  been  largely  ignored,  the  errors  Introduced  are  small  in  the 
regions  of  interest.  Devices  of  55^  efficiency  are  clearly  impractical, 
while  a  20f}  efficient  device  can  probably  be  bvillt  with  present  techniques 


and  materials. 
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Further  calcxalation  and  experiments  dealing  with  the  lifetime 
expected  from  such  devices indicate  that  operation  can  be  expected 
for  the  order  of  many  years  for  tungsten  electrodes  and  about  a  year 
for  molybdenum  electrodes. 

Data  relating  effective  work  functions  as  a  function  of  cesium 
pressure  which  are  necessary  for  calculations  of  the  type  described  are 
extremely  scarce  and  relatively  unreliable.  The  obtaining  of  more  extensive 
and  more  acctirate  information  would  allow  the  desi®i  of  thermionic  converters 
to  be  made  in  a  relatively  routine  fashion. 

Further  data  describing  the  effect  of  a  partial  cesivm  coverage 
on  "toe  thermal  emission  properties  of  surfaces  is  also  of  great  importance. 

As  has  been  already  mentioned,  the  model  of  an  ohmic  plasma 
is  crtide)  it  is,  however,  the  only  workable  model  available.  It  is  there¬ 
fore  advisable  to  derive  theoretically  and  to  verify  experimentally  a  more 
accxirate  description  of  the  effects  of  the  partially  ionized  plasma  on  the 
transport  of  electrons .  Further,  it  is  important  to  obtain  a  better  picture 
of  how  the  presstire  of  this  partially  ionized  plasma  effects  the  electronic 
space  charge  and  the  final  potential  distribution  in  the  diode. 

Causes,  control  and  uses  of  the  arcing  mode  are  little  understood. 
In  order  to  intelligently  design  a  thermionic  diode,  the  designer  must  know 
how  to  suppress  the  arc  if  it  is  not  desired  or  how  to  initiate  it  if  it 
appears  that  operation  in  that  mode  is  desirable. 
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It  is  possible  that  a  criterion  other  than  efficiency  shovild 
be  used  for  a  particular  application.  For  example,  in  a  satellite  power 
supply,  maximum  power  to  weight  ratio  would  be  requied.  Such  calculations 
can  be  done  in  a  manner  directly  analogous  to  that  described  here  if  the 
total  weight  of  the  system  can  be  calculated.  Since  the  thermionic  elements 
themselves  are  of  negligible  relative  weight,  this  means  that  about  every  oper¬ 
ating  point  of  the  diode  a  complete  power  system  must  be  designed  to  maximize 
the  power  to  weight  ratio.  This  is  at  present  an  impossible  task.  However, 
on  the  basis  of  some  simple,  not  unreasonable,  assumptions,  calculations 
have  been  carried  out  and  published^^^^  for  a  particular  system  of  a  nuclear 
reactor  and  radiator.  In  that  publication,  it  was  demonstrated  how  the  ratio 
of  output  power  to  radiator  weight  alone  could  be  maximized. 

Other  criterion  could  be  power  per  volume,  peak  power  density, 
etc.  However, regardless  of  the  criterion  used,  the  problems  pointed  out  here 
exist  and  must  be  solved  before  a  better  model  of  the  thermionic  diode 


can  be  constructed 


) 
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Explanation  of  Table  1 

Many  cf  the  figures  presented  in  the  table  (which  is  a  direct 
reproduction  of  the  computer  output)  are  in  what  is  known  as  "E- notation. " 
That  is,  a  number  written  as  3*^5E-05  =  x  lo”^  or  S.k'fE  01  =  2.47  x  10^, 
etc. 

There  are  10  columns  in  the  table,  each  corresponding  to  a  diode 
with  the  anode  temperature  given  on  the  first  line.  The  values  printed  in  the 
heading  of  the  table  pertain  to  all  10  of  the  diodes. 

Defined  in  terms  of  the  qtiantities  in  this  report,  the  symbols  in 
the  left-hand  column  are  as  follows.  Becaijse  of  the  calculations  having 
been  done  for  a  unit  area  of  cathode,  the  units  of  the  tabulated  quantities 
may  not  appear  consistent;  but  they  are: 

TEMP-A  =  T_  a  anode  temperature  (deg  K) 

Q 

TLOAD  =  temperature  of  load  which  in  these  calculations  and  in  the  paper  was 
assvimed  eqxial  to  T^^. 

CHPT-A  =  't  =  anode  work  function  (volts) 
a 

2 

SATJ-A  =  J  =  anode  saturation  current  density  (amp/cm  ) 

&s 

2 

NETSAT  =  J  -  J  =  net  saturation  current  density  (amp/cm  ) 
cs  &s 

EMIS-A  =  =  thermal  emissivity  of  anode 

VIEWFN  =  F  =  view  function 

2 

■RADHET  =  W  =  heat  radiated  thermally  $JBtts/cm  ) 

OCVLT  =  =  open  circuit  voltage  (volts) 

FLTVLT  =  4'  -  It  =  internal  voltage  in  flat  mode  (volts) 

c  a 
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INTVLT  =  V  =  internal  voltage  (volts) 

p 

FCRCUR  =  =  forward  current  density  ( amps /cm  ) 

•  p 

BAKCUR  =  Jq  =  back  current  density  ( amps /cm  ) 

p 

NETCUR  =  J  =  net  current  (amps /cm  ) 

•  2 

RUy*’  =  Rt  =  optimim  lead  resistance  for  1  cm  cathode  (ohms) 

L 

PLVDRP  =  JRp  =  voltage  drop  in  plasma  (volts) 

LDVDRP  =  JRj^  =  voltage  drop  in  optimized  leads  (volts) 

EXTVLT  =  ’^load  ~  sx'temal  voltage  applied  to  optimized  load  (volts) 

CONI^  =  =  thermal  conductivity  of  optimized  lead  (watt-cm/deg  K) 

LDTHC  =  Kj^  (Tq  -  T^)  =  heat  lost  by  thermal  conduction  down  optimized 

p 

leads  (watts /cm  ) 

PLTHJiC  =  Kp  (Tq  “  *^3)  =  lost  by  thermal  conduction  across  the  plasma 

p 

(wat ws/cm  ) 

MU  =  4)  =  (volts) 

QC  =  (d)  +  2KT^)  (watts/cm^) 

^  =  j  ($  +  2KT„)  (watts /cm^) 

COOL  =  -  Qa  -  JV  (watts/cm^) 

C  0 

JOUIEP  =  1/2  Rp  =  Joule  heating  of  plasma  which  is  retvirned  to  cathode 

p 

(watts /cm  ) 


p  ^ 

JOUIEL  =  1/2  J  R  =  Joule  heating  of  optimized  leads  which  is  returned  to  cathode 
L 

2 

(watts /cm  ) 
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POtflN  =  =  heat  in  (watts /cm^) 

p 

POWOUT  =  P  .  =  power  out  (watts /cm  ) 
0\1*C 


EFF  = 

U 

=  efficiency 

uphill  \ 

MODE  =  ( 

F 

D 

^  if  optimum  occurs  in^ 

flat  \ 

downhill  1 
<  / 

mode 
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D'Va  I94A566 


Load 


Fig.  I 

Basic  thermionic  diode 


<“)  (b)  (c) 

Uphill  Downhill  Flat 


Fig.  2 

Mode  of  operation 
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CURVE  617438 


Fig.  3 

P  as  a  function  of  y  and  the  ideal  efficiency,  tIq- 
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CURVE  517439 


Optimum  values  of  y ,  and  p  as  functions  of  the  ideal  efficiency, 


Effective  Work  Function  (volts) 
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CURVE  5.3866 


Surface  Temperature  (\) 

Fig.  7 

Work  function  as  a  function  of  cesium  bath  temperat\are  and  siurface 

temperature 


Cathode  Temperature  ("K) 
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600  700  800  900  1000  1100  1200 
Anode  Temperature  (®K) 


Fig.  9 

Efficiency  as  a  function  of  cathode  and  anode  temperatures 
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Spacing  (microns) 

Fig.  11 

Efficiency  as  a  function  of  spacing  and  cathode  temperature 
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Spacinq '  2. 5  p 


Efficiency  as  a  Function  of  Spacing  and  Anode,  Cathode,  and 
Cesium  Temperatures. 
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Fig.  12b 


Anode  Temperature  ("K ) 
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Fig.  12c 


Anode  Temperature  ("K) 
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Anode  Temperature  (®K) 


Fig.  12  e 


Anode  Temperature  ("K) 
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Anode  Temperature  (®K) 
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Fig.  12h 
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Spacing  -  30  m 


Anode  Temperature  (®K) 
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Fig.  12 J 


Anode  Temperature  (“K) 
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Fig.  12k 


Anode  Temperature  {‘’K) 
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Spacing  - 125  m 
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Anode  Temperolure  {*K ) 
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Fig.  12n 


Anode  Temperature  (®K) 


CUMVI  St04«l 


Spacing  -  25011 


Cesium  Temperature  (°K) 


Fig.  12o 
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Cathode  Temperature  (®K) 


Fig.  15 

Optimum  efficiency  as  a  function  of  cathode  temperature 


.7 


CUMVCSIftTt 


Cathode  Temperature  (®K) 


Fig.  Ik 

Optimum  anode  and  cesixsa  temperatures  as  a  function  of  cathode  temperature 
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